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ABSTRACT 


This  thesis  examines  parametric  x-radiation  (PXR)  which  is  the  Bragg  scattering  of  the 
virtual  photons  associated  with  the  Coulomb  field  of  relativistic  charged  particle  from  the 
atomic  planes  of  a  crystal.  Higher  order  parametric  x-radiation  from  the  {002}  planes  of  a 
thick  mosaic  graphite  crystal  have  been  observed.  The  raw  PXR  data  was  collected  using 
a  SiLi  detector  and  a  Pulse  Height  Analyzer  (PHA)  software  program,  fhe  data  was 
corrected  for  various  effects  including  attenuation,  detector  drift  and  efficiency  The 
absolute  number  of  photons  per  electron  was  obtained  by  using  the  fluorescent  x-ray  yield 
from  a  tin  foil  backing  on  the  graphite  crystal  to  determine  the  LIN  AC  current  The 
number  of  photons  per  electron  observed  greatly  exceeds  the  expected  values. 
Comparison  of  the  ratio  of  intensity  of  a  given  order  to  the  first  order  [I(n)/I(l)]  to  the 
theoretical  ratio  shows  that  the  ratios  increase  with  order.  Not  only  is  the  absolute 
intensity  greater  than  expected,  but  the  higher  orders  (compared  to  the  first  order)  are 
larger  than  expected.  Lastly,  the  intensity  for  various  crystal  angle  orientations  and  a  fixed 
detector  angle  was  measured. 
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I.  INTRODUCTION 


Parametric  X-radiation  (PXR)  is  the  scattering  of  the  virtual  photons  associated  with 
the  Coulomb  field  of  a  relativistic  charged  particle  by  the  atomic  planes  of  a  crystal  at  or 
near  the  Bragg  condition.  It  is  monochromatic  and  spectrally  intense  Potentially,  PXR 
may  be  an  important  new  X-ray  source  for  a  host  of  applications  [Ref  1].  Ter-Mikaelian 
(1972)  first  developed  the  theory  of  PXR  as  resonant  radiation  produced  in  a  thin  crystal 
[Ref  2],  Later  Feranchuck  and  Ivashin  developed  a  theory  for  thick  crystals.  In  1985,  the 
first  experiment  verifying  this  theory  was  performed  by  Baryshevsky  et  al  at  the  Tomsk 
synchrotron  [Ref  3].  This  experiment  used  a  900  MeV  electron  beam,  bombarding  a 
diamond  crystal  at  a  Bragg  angle  of  45  degrees  producing  x-rays  in  the  5  to  25  keV  range. 

The  first  experiments  to  observe  PXR  outside  the  USSR  were  conducted  at  the  Naval 
Postgraduate  School  (NPS).  [Ref  2].  The  first  experiments  were  done  with  20}im  Si 
[Ref  2],  44,  320  ^m  Si  [Ref  4]  and  Carbon  [Ref  5]. 

This  thesis  is  an  in  depth  study  of  the  enhancement  of  the  higher  order  PXR  due  to  the 
production  and  subsequent  attenuation  of  x-rays  in  a  thick  mosaic  carbon  graphite  crystal. 
This  work  presents  measurements  of  absolute  PXR  production  efficiency  and  the  effects 
of  experimental  apertures  on  the  bandwidth  of  PXR  spectral  line  widths,  and  makes 
comparisons  with  calculations.  Effects  due  to  beam  crystal  alignment  and  due  to  limited 
detector  aperture  field  of  view  have  also  been  investigated.  These  crystals  reveal  greater 
relative  production  rates  for  the  higher  order  x-rays  than  can  be  explained  by  target 
thickness  effects.  Furthermore,  these  measurements  show  that  PXR  production  for  all 
orders  in  a  thick  mosaic  graphite  crystal  greatly  exceed  predictions  from  calculations  based 
upon  current  theoretical  understanding. 


1 


n.  THEORETICAL  BACKGROUND 


Parametric  x-radiation  (PXR)  is  one  of  several  types  of  x-radiation  that  are  produced 
by  the  interaction  of  a  relativistic  charged  particle  propagating  in  structured  media.  These 
include  bremsstrahlung,  transition  radiation  (TR),  channeling  radiation  (CR)  and 
Che-enkov  radiation.  These  various  mechanisms  are  discussed  and  compared  in  reference 
6,  which  claims  that  PXR  is  the  most  effective  source  of  quasi-monochromatic  and 
focused  x-rays,  and  reference  7  which  compares  PXR  to  synchrotron  radiation. 

A.  PARAMETRIC  X-RAYS  (PXR) 

Parametric  x-rays  are  produced  when  the  electromagnetic  field  of  the  relativistic 
charged  particle  interacts  with  the  periodic  electric  susceptibility  associated  with  the 
crystal  lattice  structure.  When  the  charged  particle  beam  satisfies  the  Bragg  condition 
with  the  crystal  lattice,  virtual  photons  are  diffracted  as  real  x-rays.  [Ref  2] 

PXR  is  produced  when  the  photon  index  of  refraction  n(if,©)  is  greater  than  unity 
when  the  Bragg  condition 

iP={x,  +  T)\  (|^  =  |^b|)  (0 

is  fulfilled,  where  T  is  the  reciprocal  lattice  vector,  k  is  the  wave  vector  of  the  emitted 
photon,  and  O)  is  its  frequency.  This  condition  is  illustrated  in  Figure  1  As  a  result,  the 
condition  for  Vavilov-Cherenkov  radiation 

1-— n(ic.(a)  =  0,  (ft  =  c  =  l)  (2) 

w 

is  satisfied  in  a  crystal.  This  is  in  contrast  to  a  homogeneous  medium,  where  the  index  of 
refraction  is  less  than  unity  for  x-rays,  and  Vavilov-Cherenkov  radiation  is  not  produced. 
Therefore,  PXR  is  simply  Vavilov-Cherenkov  radiation  produced  by  diffraction  according 
to  equation  (1).  [Ref  8] 
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(reciprocal  lattice  vector) 

Crystal  planes 


parametric  x-rays  electromagnetic  field  of 

(real  photons)  charged  particle  beam 

(virtual  photons) 

Figure  1 .  Diffraction  of  virtual  photons  associated  with  the  charged  particle  beam  by  the 
atomic  planes  of  the  crystal  showing  the  Bragg  angle  as  the  angle  between  the  charged 
particle  beam  and  the  crystal  planes.  [Ref  2] 

There  are  two  ways  to  treat  the  diffraction  of  real  and  virtual  x-rays:  1 )  dynamic  theory 
and  2 1  kinematic  theory.  The  first  case  is  realized  in  ideal  crystals  of  thickness  L,  greater 
than  the  X-ray  extinction  length.  On  the  other  hand,  the  kinematic  theory  of  PXR  is 
applicable  for  real  crystals,  which  consist  of  thin  mosaic  blocks  turned  relative  to  each 
other  at  an  angle  8>tn/E,  where  m  is  the  mass  of  the  charged  particle  and  E  is  its  energy. 
The  kinematic  diffraction  is  more  convenient  for  the  experimental  observation  of  PXR  and 
it  can  be  most  easily  analyzed  because  the  angular  and  spectral  distributions  of  PXR  are 
essentially  simplified  and  have  a  universal  form  for  different  crystals.  [Ref  8] 

Ter-Mikaelian  formulated  the  expression  for  the  number  of  photons  emitted  by  a 
charged  particle  in  a  thin  crystal;  treating  the  interaction  as  a  perturbation,  with  the  crystal 
requirement 

>cL|n-lj<l,  (3) 

where  L  is  the  effective  thickness  of  the  crystal.  If  the  refraction  of  photons  and  the 
multiple  scattering  of  the  electrons  in  the  crystal  are  taken  into  account,  then  the  results. 
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for  a  mosaic  crystal  with  8>ni/E  and  with  PXR  from  different  blocks  being  coherent, 
remain  correct.  The  spectral  and  angular  distributions  of  emitted  photons  are  formulated 
by 


— —  =  — oi,.  gJ 
dn^dci)  In  * 


Ik’j^.CUV  +  tj" 


_  I -exp  5((?). 

\  / 


^  =  -*fr  +  r,-y[Re{go)-0,].  (5) 

Here  g-j  and  go  are  the  Fourier  components  of  the  dielectric  susceptibility  parallel  and 

perpendicular  to  the  reciprocal  lattice  vector.  These  values  are  directly  related  to  the 
coherent  scattering  amplitude  of  the  photons  by  the  atoms  of  the  crystal  Bg  is  the  angle  of 
multiple  scattering.  Equation  (4)  is  written  in  a  coordinate  system  with  the  z-axis  directed 
along  the  particle  velocity  v;  =  con^  ,  and  L,  =:[coIm{g5)]  [Ref  8]  The  differential 

radiation  cross  section  for  PXR  as  derived  by  Feranchuck  and  Ivashin  can  also  be  written 
in  terms  of  angular  displacements  from  the  Bragg  angle  instead  of  as: 


1  eM  ,  ,2  ^,cosM20j  +  0: 


dd^dd^do)  he  (Og  +  ^  +  +  g, 


16sin  Bg  — - +p 

V  tanBj  Q)g 


where  g^  =  n(to)-l=—  —  ,  n(a))  is  the  index  of  refraction,  (Op  is  the  plasma 

21  oj  J  Z  ^ 


frequency,  = - —  ,  where  angular  displacements  away  from  the 
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Bragg  angle  in  and  out  of  the  incident  plane  respectively,  |gj  =  ---- -  /,  j  fj  js  the 

x-ray  scattering  factor  of  the  crystal,  V  is  the  volume  of  the  unit  cell  of  the  crystal,  tg  is 
the  classical  electron  radius,  ^i=c=l,  p^=Xl{L^{oi)2it)  and  La(o))  is  the  photon 

absorption  length,  [Ref  8] 

The  characteristic  quasi  monochromatic  photons  produced  at  large  angles  relative  to 
the  charged  particle's  velocity  vector,  v,  distinguishes  PXR  from  other  radiation 
mechanisms.  The  particle  energy  does  not  determine  the  emitted  photon's  wave  vector,  k, 
nor  frequency  (03.  The  emitted  photon's  wave  vector  depends  upon  the  reciprocal  lattice 
vector,  T,  the  velocity  vector  of  the  electron  and  the  Bragg  frequency  The  Bragg  angle, 
is  the  angle  between  the  velocity  vector,  v  and  ^he  crystallographic  planes  associated 
with  the  reciprocal  lattice  vector  r.  The  PXR  maxima  are  defined  by  the  crystal  lattice 
structure  amplitudes  corresponding  to  T  The  values  of  the  emitted  wave  vectors,  ic,  and 
frequencies  0)3  are 

k=0)gV*Z;  jc’v  =  2Bg,  (7) 


a)-(n)  = - =  — — —  = - ,  n=  1,2,3  .. 

2|t,|  2sin6j  dsindg 


where  d  is  the  distance  between  the  crystal  planes  corresponding  to  the  vector  T  [Ref  3  .] 
The  intensity  of  PXR  is  proportional  to  jg^f  and  is  greater  in  crystals  with  higher 


packing  factors,  such  as  a  diamond  lattice.  Rotation  of  a  crystal  by  some  A0  results  in  the 
emitted  wave  vector,  k,  rotating  2A0  in  the  same  direction  while  the  frequency  varies  as  a 


result  of  equation  (7),  The  angular  distribution  of  emitted  photons  can  be  derived  by 
integrating  equations  (4)  or  (6)  over  o),  resulting  in 


d^N  _  y  CJfl  ^  ^  cos- 20^  +  0;] 

de,de,  ^,4rt  c  \  ^aj  sin‘0fl 
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(10) 


/\2  ^  /-\2  I  I 

^ph  =p-+®i 
^  =  Cultv=a.  +  Cam-d.v  ■*■ 


(11) 


where  m  is  the  rest  mass  of  the  electron,  is  the  absorption  coefficient  which  is  a 
function  of  o)  and  0ph  is  the  characteristic  angular  spread  of  the  PXR.  0s  includes  the 
effects  of  multiple  scattering  of  the  beam  electrons,  electron  beam  divergence  and  crystal 
mosaicity.  This  can  be  put  in  dimensionless  form  by  means  of  the  normalized  amplitude  J= 
N/Nq  and  angles  x,y  s  0x,y(0ph> 


dxdv 


,/  V  jr'cos' 20.  +  V' 
(jc'+y-  +  lj 


i-exp - 


sin  0. 


(12) 


(13) 


The  number  of  photons  detected  in  an  angular  width  of  0^  about  0^  is  obtained  by 
integrating  equation  9: 


Nc  =  7tN,{l  +  cos-  20a  jJJ'’ 


A/’o  =  iV,{l  +  cos^  208)|ln- 


(14) 


(15) 


0 

where  .  The  detector  angular  size  is  a  factor  in  determining  the  value  of  Nq, 

even  for  0D*'^®ph  factor  does  not  occur  in  channeling  or  in  bremsstrahlung 

radiation;  however  in  PXR,  the  intensity  slowly  decreases  inversely  proportional  to  0  as  0 
increases.  [Ref  8] 
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According  to  equations  (14)  and  (15),  there  is  no  threshold  characteristics  in  the  PXR 
intensity  The  angular  frequency,  however,  depends  on  6ph(E)  which  increases  with  the 
electron  energy 


E<E  , 

Opt 


I  m '  +  0; 

V  Isoi  • 


(16) 


in  accordance  with  equation  (9)  and  (10),  thus  setting  the  limiting  threshold  energy  for 
PXR.  The  number  of  photons  decreases  rapidly  for  E<EQpt  The  number  detected 
becomes 


N„s|w,(l  +  cos-29,)i 


fi^l 

4 

^  B 

^  ®of>t  ^ 

(H) 


This  phenomena  suggests  that  if  the  beam  energy  is  not  sufficient,  regardless  of  whether 
all  the  other  parameter  conditions  are  met,  PXR  may  not  be  observed.  [Ref  8] 

Integration  of  equations  (4)  or  (6)  over  or  dQ  respectively  results  in  the  frequency 
distribution; 


(18) 


u  =  tanfrg 


0)b^p>, 


(19) 


Equation  ( 1 8)  can  be  put  in  terms  of  frequency  by  usin®  the  chain  rule 


dN  _  dN  du 
doi  du  dco 


and 


—  can  be  determined  using  equation  (19).  Although  equations  (18)  and  (19)  do  not 
do) 

apply  to  real  experiments  with  finite  aperture,  they  have  been  previously  applied  to  the 
experimental  results.  [Ref.  8] 
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B.  MULTIPLE  SCATTERING,  BEAM  DIVERGENCE  AND  MOSAICITY 


In  reference  1,  the  Adischev  et  al.  have  used 


e 


sts 


=  L  |_L 


(20) 


Es=21  MeV,  L  is  the  effective  thickness  of  the  crystal  and  Lr  is  the  radiation  length  This 
value  is  an  approximation  applicable  to  very  high  energy  electrons  and  is  well  known  to 
over  estimate  the  multiple  scattering  effect.  [Ref  9]  The  Bethe-Ashkin  formula  for 
multiple  scattering  which  gives  a  more  accurate  estimate  for  0^  is  used  in  the  analysis 


and  is  given  by; 


A  {pv)‘ 


(21) 


where  Z  is  the  atomic  number,  A  is  the  atomic  weight  and  t  is  the  thickness  of  the  foil  in 
grams  per  square  centimeter.  The  charge  of  the  scattered  particle  is  z,  p  is  its  momentum 
and  V  its  velocity,  <02>  is  an  average  scattering  angle  value  obtained  by  integrating  over 
angle,  where  0i  is  the  maximum  angle  used  in  the  analysis.  It  is  estimated  as  the  angle  for 
which  there  is,  on  the  average  throughout  the  full  thickness  t  of  the  foil,  only  one  collision 
with  0  >  01  [Ref  11] 

Since  x-rays  are  produced  throughout  the  electron's  path  length,  the  amount  of  multiple 
scattering  varies  from  zero  to  the  value  appropriate  for  the  full  path  length  of  the  electron 
in  the  crystal  Additionally,  the  x-rays  which  are  generated  inside  the  crystal  are 
attenuated  due  to  absorption.  Thus  there  is  a  weighting  of  x-rays  reaching  the  detector, 
with  the  weighting  factor  equal  to  «■*“,  where  ^  is  the  absorption  coefficient  for  the 
particular  x-ray  energy.  Therefore,  x-rays  produced  near  the  front  face  of  tiie  crystal  are 
weighted  more  than  those  near  the  back  surface.  This  weighting  effectively  biases  the 
detected  x-rays  toward  those  produced  at  lower  values  of  x,  i  e  near  the  front  face  of  the 
crystal  where  the  multiple  scattering  is  lowest.  The  proper  way  to  account  for  the  effects 
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of  multiple  scattering  would  be  to  add  up  all  the  x-rays  produced  as  a  function  of  x, 
including  an  x  dependent  multiple  scattering  effect. 

For  simplicity,  a  mean  value  of  can  be  calculated  by  weighting  the  value  of 

<02>  with  the  absorption  attenuation  of  x-rays,  The  mean  value  for  x  is  defined  as 
follows 


(22) 


[Ref  9]  Since  the  Bethe-Ashkin  formula  for  <0j^  §  2>  iq  first  order  is  proportional  to  x, 
this  value  for  x  is  used  in  equation  (21)  in  place  of  the  crystal  thickness,  t,  to  determine  the 
multiple  scattering  angle  0s  [Ref  11].  In  our  treatment,  the  effect  of  the  divergence  of 
the  electron  beam  is  not  included  To  more  rigorously  include  the  effect  of  beam 
divergence  on  the  x-ray  angular  distribution  and  bandwidth,  the  electron  beam  angular 
distribution  must  be  convolved  with  the  an^lar  distribution  of  the  x-ray  production  as  has 
been  done  e  g.  in  the  case  of  transition  radiation  [Ref  12],  For  simplicity  of  calculation 
the  effect  of  divergence  was  not  included  in  the  analysis  since  the  effective  beam 
divergence  of  the  NFS  LINAC  (<1  mrad)  is  unimportant. 


C  EFFECTS  OF  FINITE  APERTURE 

Previous  PXR  experiments  to  date  have  compared  the  experimental  data  with  the 
results  of  equation  (18).  However,  the  theoretical  values  obtained  by  using  equation  (18) 
refer  to  an  unrealistic  idealization  i.e.,  an  infinite  detector  aperture,  which  does  not 
approximate  the  real  conditions  of  these  experiments  and  neglect  the  effect  of  the  aperture 
on  the  PXR  spectrum. 

In  order  to  calculate  the  spectral  distribution  for  real  experimental  conditions,  equation 
(15)  should  be  integrated  over  the  appropriate  angular  field  of  view  as  determined  by  the 
detector  size  and  the  distance  from  the  crystal  to  the  detector.  The  integration  limits  6, 
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and  6^  are  set  by  the  angular  fields  of  view  in  and  out  of  the  plane  of  observation 
(A0^  =  20,  and  A0,  =20,  respectively).  [Ref  13] 

Case  I  Infinite  0,,  finite  0,  Integrating  equation  (7)  first  over  0y  with  a  finite  upper 
Untit  of  0,  and  then  integrating  0^  from  0  to  <»  results  in 
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where  the  function  J2  is  defined  as; 
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and  ttj,  = 


This  is  derived  by  using  a  delta  function  approximation  for  the  x 


dependence,  see  equation  6.  [Ref  9] 

Case  11  0,  and  0^  both  finite.  In  this  case  there  are  two  regimes,  a)  The  first  is  the 

narrow  line  limit  regime  in  which  p,  =  A0,,,  then 


d(0 


(25) 


b)  The  second  is  the  aperture  limited  bandwidth  regime  for  which  the  delta  function 
approximation  is  still  correct  if  A0„  »  p.  Then: 


In  this  case  the  energy  bandwidth  is  determined  based  on  the  relative  value  of  A6^  and 
the  full  width  at  half  maximum  (FWHM)  of  If  A0,  is  larger  than  the  FWHM  of 

then  the  bandwidth  is  determined  to  be  the  FWHM  of  J^iu).  This  region  is 
referred  to  later  in  this  thesis  as  the  "near  region"  If  Ad^  is  smaller  than  the  FWHM  of 
y,(«),  then  the  bandwidth  is  proportional  to  Ad^  and  is  relatively  insensitive  to  the  form 
of  Jiiu)  This  regime  is  called  the  small  aperture  regime  referred  to  later  in  this  thesis  as 

the  "far  region".  The  bandwidth  for  the  small  aperture  regime  is: 


Aa)  = 


tan{0e) 


It  is  the  small  aperture  regime  that  applies  to  this  thesis  work  and  that  is  most  applicable  to 
previous  PXR  experiments  [Ref  10],  Figures  2  and  3  illustrate  the  function  J^iu)  with 

A0^<FWHM  of  y3(M)for  the  far  region  and  Ad_,>FWHM  of  J.iu)  for  the  near  region 
respectively. 
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Figure  3 .  J2  and  the  field  of  view  for  the  "near  region”  at  the  Bragg  angle. 
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ra.  PXR  EXPERIMENTAL  DESCRIPTION 


Several  PXR  experiments  were  conducted  using  a  I  .39±0.0i  mm  thick  graphite  crystal 
in  Bragg  geometry.  The  x-ray  observation  angle  was  set  at  45  degrees  with  respect  to  the 
electron  beam  direction.  In  this  configuration,  corresponding  to  a  Bragg  angle  of  22.5 
degrees  with  respect  to  the  {002}  atomic  planes  of  the  carbon  mosaic  crystal, 
measurements  were  made  at  detectors  placed  at  100  cm.  and  29  cm.  from  the  target  to 
observe  the  effects  of  changing  the  solid  angle  of  observation.  At  both  these 
configurations,  "rocking  data"  measurements  were  made  "Rocking  data"  were  obtained 
by  rotating  the  crystal  about  the  nominal  Bragg  angle  of  22.5  degrees.  The  rocking  data 
provided  the  dependence  of  the  absolute  and  relative  intensities  of  the  various  order  x-rays 
when  the  Bragg  condition  was  not  exactly  satisfied. 

A.  ALIGNMENT  AND  ELECTRON  BEAM  STEERING 
The  NPS  Linear  Accelerator  (LINAC)  was  used  to  conduct  these  experiments.  The 
beam  energy  varied  with  each  experiment  with  a  nominal  electron  energy  of  90  MeV. 
Two  different  SiLi  detectors  [Ref  14  &  15]  were  used  in  these  experiments.  The  first  was 
an  EG&G  ORTEC  model  7113-16250  with  an  active  detector  region  of  200  mm^.  and 
the  other  was  a  Canberra  model  SI200250  with  an  active  detector  region  of  200  mm^. 
The  SiLi  detectors  have  nominal  time  resolution  times  of  12  nsec  determined  by  the  pulse 
shaping  preamplifier.  Because  the  LINAC  is  pulsed  at  60  Hz  with  a  macro  structure 
length  of  1  |isec,  the  system  is  count  rate  limited.  The  LINAC  duty  cycle  makes  the  count 
rate  60  MHz  when  a  single  photon  is  detected  during  each  LINAC  pulse.  The  LINAC 
was  adjusted  to  limit  the  macroscopic  average  count  rate  including  background  to  one 
count  per  three  to  five  machine  macropuises  in  the  whole  detector  energy  range. 
Consequently,  the  LINAC  beam  was  limited  to  a  dark  current  estimated  to  be  of  the  order 
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of  2x10"  amps.  Dark  current  refers  to  using  the  LIN  AC  with  zero  gun  grid  voltage,  so 
that  only  stray  electrons  are  accelerated  The  ratio  of  counts  per  machine  pulse  was 
determined  by  use  of  fluorescent  foils  to  prevent  the  double  counting  of  photons  by  the 
detector  That  is,  if  two  photons  are  stopped  in  the  detector  simultaneously  (within  the 
peak  shaping  time  of  the  detector)  these  two  photons  will  appear  as  a  single  photon  with 
an  energy  equivalent  to  the  sum  of  the  two  coincident  photons.  Pulse  pile  up  was  of 
critical  concern  in  this  experiment  because  two  first  order  PXR  photons  observed 
simultaneously  would  register  as  a  second  order  PXR  photons  at  double  the  energy  This 
would  falsely  enhance  the  higher  orders  as  well  as  reduce  the  first  order,  resulting  in  ratios 
of  intensity  of  higher  order  peaks  to  first  order  peaks  too  high  Maintaining  a  constant 
dark  current  was  difficult  to  achieve  and  required  constant  operator  attention. 

The  experimental  setup  is  shown  in  Figure  4  The  beam  is  shown  entering  a  vacuum 
chamber  from  the  left  at  the  entrance  port.  At  the  center  of  the  chamber  is  a  target  ladder 
which  can  be  remotely  rotated,  raised,  lowered  and  tilted.  An  alignment-reference  laser  is 
located  outside  the  port  located  90  degrees  counter  clockwise  from  the  entrance  port 
The  detector  is  located  45  degrees  clockwise  from  the  beam  exit  port  or  135  degrees 
counter  clockwise  from  the  entrance  port.  To  help  align  the  electron  beam  a 
phosphorescent  screen  was  placed  at  the  end  of  a  27.5  in.  pipe  located  at  the  beam  exit 
port. 

Experiments  prior  to  July  22  1992  were  conducted  using  an  EG&G  ORTEC  Si(Li)  x- 
ray  detector  model  number  7113-16250.  Due  tc  vhe  fixed  height  of  the  target  chamber 
and  the  detector's  large  liquid  nitrogen  tank  located  below  the  detector  window  the 
detector  couldn't  be  placed  closer  than  about  two  feet  from  the  chamber.  For  experiments 
with  this  detector  a  copper  vacuum  pipe  27.5  in  long  was  placed  at  this  port  fitted  with  a 
one  mil  Kapton  window  (1  in.  diameter).  The  results  obtained  frcm  these  experiments 
will  be  referred  to  as  the  far  case  or  the  far  region.  Experiments  on  or  after  July  22'^<^ 
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Figure  4.  PXR  experimental  setup  in  the  LINAC  end  station,  showing  incoming  charged 
particle  beam  and  beam  line,  alignment-reference  laser,  detector,  and  entrance  and  exit 
ports.  1.  Target  ladder  with  calibration  foil,  mirror,  and  carbon  crystal.  2.  downstream 
phosphorescent  screen. 

were  conducted  using  a  Canberra  Si(Li)  x-ray  detector  model  number  SI200250,  This 
detector  consist  of  a  two  fill  and  vent  port  dewar  to  allow  operation  at  any  orientation  and 
a  "slim-line”  detector  chamber  which  is  small  in  size  and  weight.  [Ref  14]  This  detector 
has  a  much  smaller  horizontally  mounted  liquid  nitrogen  tank  and  could  be  positioned 
directly  next  to  the  vacuum  chamber.  The  vacuum  chamber  port  was  fitted  with  the  same 
one  mil  Kapton  window  that  covered  the  copper  pipe.  The  results  obtained  from  these 
experiments  will  be  referred  to  as  the  near  case  or  near  region. 

A  remotely  vertically  positionable  ladder  and  2  axis  goniometer  were  used  to  hold  and 
orient  the  mosaic  carbon  crystal  used  in  this  experiment.  The  crystal  used  in  ail  the 
experiments  is  a  mosaic  graphite  monochromator  measured  to  be  1.39  ±  0.01  mm.  thick, 
manufactured  by  Union  Carbide.  The  density  of  the  crystal  is  2.260  ±  0.005  grams  per 
cubic  centimeter.  The  spacing  of  the  reflecting  (002}  planes  is  3.356  ±  0.003  angstroms. 
The  mosaic  spread  is  ^.4  degrees.  (Ref  16]  (Pure  crystals  (mosaic  spread  =  0)  such  as 
Si  have  lower  reflectivity  than  mosaic  crystals  in  normal  x-ray  diffraction  and  therefore,  a 
higher  the  mosaic  spread  results  in  a  theoretically  higher  measured  intensity.)  On  June  5* 
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half  of  the  crystal  was  covered  with  a  0.0275  mm  thick  layer  of  tin  7.3 1  g/cm^  on  the  back 

face  of  the  crystal.  This  foil  was  essential  as  a  calibration  source.  From  the  location  of  the 
fluorescence  line,  the  energy  stability  and  energy  calibration  of  the  detector  through 

out  the  experiment  could  be  monitored.  The  observed  yield  of  photons  is  proportional  to 
the  integrated  electron  current  through  the  foil.  This  information  was  used  to  determine 
the  absolute  photon  production  rate.  [Ref  1 7] 

The  ladder  which  housed  the  mosaic  graphite  crystal  also  contained  three  separate  foils, 
copper  (Cu),  yttrium  (Y)  and  indium  (In),  which  were  replaced  on  June  5^  by  a  sandwich 
foil  consisting  of  a  titanium  (Ti)  foil  0.0164  mm.  thick,  a  yttrium  (Y)  foil  0  0603  mm. 
thick,  and  a  tin  (Sn)  foil  0.0275  mm.  thick.  These  foils  were  used  for  energy  calibration 
purposes  and  to  determine  the  LINAC  beam  current.  (The  details  of  this  calculation  are 
presented  in  reference  17  )  The  ladder  also  contained  a  phosphorescent  screen  with  a  pin 
hole  in  the  center  which  was  used  to  position  the  electron  beam.  The  target  ladder  also 
contained  a  mirror,  coplanar  with  the  crystal  face  which  was  used  to  reflect  a  reference 
laser  back  onto  itself  to  establish  the  zero  degrees  or  "home"  position. 

The  reference  laser  was  used  to  establish  the  orientation  and  position  of  the  crystal 
with  respect  to  the  electron  beam.  Using  cutouts  with  a  pinhole  in  the  center  to  cover  the 
ports,  the  laser  was  adjusted  so  that  is  went  through  the  geometric  center  of  the  chamber. 
The  ladder  was  then  placed  in  the  center  of  the  chamber  with  the  motor  controls  attached 
and  the  pinhole  set  to  the  position  of  the  laser  beam.  The  ladder  was  manually  adjusted 
for  tilt  so  that  the  laser  went  through  the  pinhole  in  the  center  of  the  phosphorescent 
screen,  and  remained  on  the  pinhole  upon  360  degree  rotation  of  the  crystal.  A  mirror  and 
the  graphite  were  positioned  so  that  their  faces  were  coplanar  and  such  that  the  rotation 
axis  of  the  ladder  was  in  their  planes  as  well.  Then,  by  setting  the  mirror  at  45  degrees 
with  respect  to  the  laser  beam,  the  laser  beam  was  reflected  down  the  exit  pipe.  A  second 
phosphorescent  screen  was  placed  at  the  end  of  the  exit  pipe  with  the  center  of  the  screen 
location  based  on  the  location  of  the  reflected  laser.  This  second  screen  greatly  improved 
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the  alignment  of  the  experiment.  It  was  determined  that  the  position  of  the  center  of  the 
screen  was  not  identically  at  the  geometric  center  of  the  pipe.  After  the  first  data  were 
obtained,  using  this  screen,  the  screen  slipped,  and  could  not  be  replace  exactly  without 
breaking  vacuum  and  re-configuring  it.  The  screen  was  replaced  with  its  center  near  the 
geometric  center  of  the  pipe,  which  introduced  a  slight  additional  uncertainty  in  the  crystal 
angle  of  less  than  two  tenths  of  a  degree.  For  this  reason,  in  future  experiments,  OTR 
should  be  used  on  a  down  stream  foil  as  an  alignment  as  was  first  done  with  Si.  [Ref  4] 

After  the  chamber  was  sealed  and  a  vacuum  was  established  the  LINAC  was  then 
tuned  to  the  given  energy,  and  a  beam  was  established.  The  beam  was  adjusted  so  that  it 
went  through  the  center  of  the  quadrupole  magnets,  so  that  adjustments  in  klystron 
frequency  and  phase  during  the  experiment  would  not  affect  the  location  of  the  beam. 
Using  end  station  steering  magnets,  the  beam  was  steered  so  that  it  went  through  the  pin 
hole  located  on  the  screen  on  the  ladder.  Since  the  only  means  of  steering  without 
changing  the  position  of  the  beam  through  the  quadrupole  magnets  is  end  station  steering 
(see  Figure  5)  a  compromise  had  to  be  made  to  get  the  beam  through  both  the  center  of 
the  screen  on  the  ladder  and  through  the  center  of  the  downstream  screen.  Additionally, 
the  beams  focus  required  a  compromise.  Focusing  downstream  leads  to  a  more  parallel 
beam  at  the  target,  but  it  also  makes  a  larger  beam  spot.  Focusing  on  the  target  gives  a 
smaller  beam  spot,  but  the  beam  electrons  are  not  parallel.  The  beam  was  focused 
approximately  half  way  between  the  two  screens  based  on  these  two  factors  This 
uncertainty  however,  contributes  to  the  uncertainty  in  the  crystal  angle  with  respect  to  the 
beam  line  of  less  than  one  tenth  of  a  degree. 

Once  the  beam  was  focused  and  the  steering  adjusted,  the  electron  gun  was  turned  off 
and  the  current  was  reduced  to  dark  current.  The  ladder's  angle  was  adjusted  using  the 
laser  since  the  elearon  beam's  position  is  known  relative  to  the  laser.  Refer  to  Figure  6 

The  angle  0c,  for  positioning  the  crystal  is  determined  with  respect  to  the  reference 
laser.  0c  is  the  angle  by  which  the  crystal  was  rotated  and  was  equal  to  0^  or  03  +  1 80 
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Figure  5.  Naval  Postgraduate  School  100  MeV  Linear  Accelerator  (LIN AC)  showing 
klystrons,  quadrupole  magnets  and  end  station  steering  magnets. 
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Figure  6.  The  relationship  between  the  incoming  particle  beam  and  the  ladder,  Og,  and 
between  the  beam  and  the  detector  which  is  fixed  at  45  degrees. 

degrees.  The  home  or  zero  degrees  position  was  established  by  using  a  mirror  on  the 

ladder.  When  the  ladder  is  oriented  at  zero  degrees  with  respect  to  the  laser,  the  beam  is 

reflected  back  onto  the  laser.  Using  a  remote  camera  the  "home"  or  zero  position  is 

established  and  fine  adjusted  after  e^'ery  rotation  of  the  crystal.  This  gives  a  more 

consistent  angle  and  also  removes  the  backlash  inherent  in  the  motor  controls.  With 
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reference  to  the  laser  the  electron  beam  enters  the  chamber  at  90  degrees  The  target 
crystal  orientation  for  these  experiments  varied  with  respect  to  the  beam  line  while  the 
detector  angle  was  fixed  at  45  degrees. 

B.  ELECTRONICS  FOR  SPECTROSCOPY 

Because  the  end  station  is  relatively  exposed  to  the  klystron  gallery  (see  Figure  5)  the 
SiLi  detectors  are  susceptible  to  ground  loops  and  radiated  noise  from  the  klystrons 
Consequently  great  effort  was  taken  to  insure  sufficient  grounding  and  shielding  to 
suppress  the  noise  pulse.  In  addition  a  gating  system  was  set  up.  Figure  7  is  a  block 
diagram  of  the  electronics  setup.  The  signal  is  received  by  either  SiLi  detector  It  is 
amplified  by  an  ORTEC  Amp  571  pre-amplifier  and  an  ORTEC  450  amplifier  which  was 
replaced  on  May  8^  by  a  TENNELEC  TC  244  amplifier.  It  is  passed  though  an  ORTEC 
426  linear  gate  with  a  gate  width  of  5  microseconds.  This  gate  was  replaced  in  later 
experiments  with  a  TENNELEC  TC  304  linear  gate  which  has  wider  gating  capabilities 
The  gate  is  triggered  by  a  Stanford  Research  System  Inc.  model  DG535  four  channel 
digital  delay/pulse  generator.  The  detector  and  preamplifier  were  located  in  the  LINAC 
end  station.  To  allow  for  possible  changes  in  settings  the  remaining  components  were 
located  in  the  control  room.  This  resulted  initially  in  a  noise  problem  that  was  eliminated 
by  trial  and  error,  using  grounding  straps  and  a  return  signal  from  the  amplifier  in  the 
control  room  to  an  oscilloscope  in  the  end  station.  The  delay  generator  delay  time  (Tj)  is 
r^justed  with  respect  to  Tq,  the  machine's  start  sequence  time.  The  delay  time  is  adjusted 
so  that  the  start  of  the  gate  coincides  approximately  with  the  apparent  arrival  of  the  beam 
pulse  at  the  ladder  in  the  chamber.  The  ideal  delay  time,  T(j  for  the  gating  was  determined 
by  using  an  oscilloscope  along  with  the  pulse  height  analyzer  and  the  fluorescent  signal 
from  a  copper  foil.  T^j  was  determined  to  be  the  time  at  which  the  maximum  signal  was 
observed.  Shorter  and  longer  delay  times  gave  no  signal.  The  maximum  was  found  to  be 


19 


X 


KliiH! 

listgj'airgti'rawgi 

■ 

■ 

Stanford  Research  Svs  Inc. 

DG  535 
pulse  generator 

I  l.PSlAG  t'  '  *" . 1  Higital  mnnn»r  I 

Figure  7.  Electronic  setup 

6  ^is,  for  the  earlier  experiments  using  the  ORTEC  amplifier,  and  30  6  |is  for  the 
TENNELEC  amplifier. 

The  LINAC  generates  60  pulses  of  one  microsecond  duration  every  second  Since  x 
rays  are  only  produced  during  the  actual  pulse  time  this  gating  doesn’t  affect  the  PXR 
spectrum  or  the  sandwich  foil  calibration  spectrum,  other  than  by  limiting  the  noise  that  is 
non-coincident  with  PXR  production.  The  gating  limits  detection  to  the  on-time  of  the 
LINAC  pulse  However  klystron  RF  noise  adds  a  bias  to  the  pulse  signal  energy 
calibration  which  tends  to  be  non-random  with  respect  to  Tq  hence  shifting  the  energy 
calibration  compared  to  that  obtained  using  radioactive  sources  Since  the  PHA 
determines  the  energy  (or  charmel)  based  on  the  energy  deposited  in  the  active  detector 
region  at  a  given  time  and  if  the  x  ray  arrives  during  a  consistently  negative  klystron  noise 
the  PHA  determines  the  energy  of  the  x  ray  to  be  lower  than  it  is.  The  klystron  noise  was 
typically  in  the  negative  voltage  portion  during  the  gating  (see  Figure  8)  The  negative 
voltage  adds  to  the  positive  voltage  produced  by  the  detected  x-rays  This  lowers  the 
output  voltage.  The  net  effect  is  a  lower  pulse  height,  which  is  interpreted  by  the  pulse 
height  analyzer  as  a  lower  energy  photon.  Therefore  the  calibration  must  be  obtained 
during  the  same  portion  of  the  machine  cycle  as  the  PXR  spectrum  to  get  the  same 
klystron  noise  effect  for  both  PXR  and  calibration  x  rays.  For  radioactive  sources  gating 
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Figure  8  Typical  shape  of  klystron  noise  during  PXR  data  collection  The  arrow 
indicates  approximate  gating  portion 

effectively  eliminates  the  signal.  Without  gating  the  noise  contributes  randomly  and 
therefore  results  in  a  net  zero  shift  Gated  fluorescent  source's  peaks  and  the  gated  PXR 
spectral  peaks  will  consistently  shift  to  lower  values  than  expected  if  the  calibration  is 
based  on  radioactive  sources  Therefore,  the  known  fluorescent  energy  values  were  used 
to  calibrate  the  detectors 

C.  ENERGY  CALIBRATION 

For  each  experiment  an  energy  calibration  spectrum  was  obtained  These  were  obtained 
by  measuring  the  x-ray  fluorescence  from  either  three  separate  foils  or  the  sandwich  foil 
described  in  Section  A,  Observed  x-rays  are  the  fluorescent  x-rays  at  the  K-edge  of  the 
excited  atom.  The  K-edge  x-ray  energy  for  titanium  is  4,509  keV,  8  04  keV  for  copper. 
14  932  keV  for  yttrium,  24,10  keV  for  indium  and  25  196  keV  for  tin  [Ref  18]  The 
calibration  is  obtained  by  establishing  the  centroid  of  the  fluorescent  peak  That  centroid 
channel  is  then  set  equal  to  the  known  de-excitation  energy  Assuming  a  linear 
relationship  between  channel  position  and  energy,  the  equation  for  determining  the  energy 
for  a  given  channel  is  determined  by  the  following  equations: 


E  =  a(chn)+/3 , 

(28) 

E,  -E, 

(29) 

a  =  — * - = —  , 

chn,  -chn^ 

/3=Ej  -a(chnj). 

(30) 

where  Ej  and  E2  are  the  energies  of  peak  I  and  peak  2  respectively,  and  chn|  and  chn2 
are  the  corresponding  centroid  channel  numbers,  a  is  the  slope  and  P  is  the  intercept  of 
the  Energy  vs,  Channel  line. 
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IV.  DATA  AND  ANALYSIS 


A.  SINGLE  ANGLE  EXPERIMENTS  (FAR  REGION, ^0;^  <  FWHM  OF  J2) 

To  obtain  a  PXR  spectrum,  the  target  ladder  was  lowered  until  the  mirror  was  in  line 
with  the  electron  beam,  and  the  home  position  was  determined  as  described  in  Chapter  III 
section  A  The  ladder  was  lowered  to  the  crystal  and  then  rotated  to  the  desired  angle 
All  single  angle  experiments  used  the  ORTEC  SiLi  detector  and  were  in  the  far  region 
For  the  far  region  the  detector  was  approximately  100  cm  from  the  target  and  subtended 
2.0x10''^  sr.  A  summary  of  the  various  variables  for  these  experiments  can  be  found  in 
Table  I,  and  the  details  are  given  below 


TABLE  I  SUMMARY  OF  VARIABLES  FOR  THE  SINGLE  ANGLE  EXPERIMENTS 
WITH  NOMINAL  ANGLES  NEAR  THE  BRAGG  ANGLE  OF  22  5  DEGREES 


Date 

Beam 

Energy 

(MeV) 

Amplifier 

(peaking  time) 

Gate 

(width) 

[delay] 

Detector 
(distance  to 
detector) 

Nominal 

Angle 

(degrees) 

5  May 

~90 

Ortcc  450 

(12  ^ls) 

Ortec  426 
(5  ms) 
[30.6  ms) 

EG&G  Ortec 
SiLi  model 
#7113-16250 
(1006  mm) 

22.5 

8  May 

~90 

*» 

If 

(1003  mm) 

22.5 

7  July 

88 

ff 

22.5 

88 

ft 

If 

23 

8  July 

84.5 

TENNELEC 
TC  244 
(8  ms) 

Ortec  426 
(35  ms) 
17.7  MS) 

ft 

(1003  mm) 

22.5 

"Beam  energy"  is  the  energy  of  the  LINAC  electron  beam,  "Amplifier"  is  the  model  of  the 
amplifier,  with  the  amplifier  peaking  time  listed  below  in  parentheses,  "Gate"  is  the  model 
of  the  linear  gate  used  with  the  gate  width  listed  in  parentheses  and  the  delay  time  from 
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the  LFNAC  start  for  the  start  of  the  gate  are  listed  in  brackets,  "Detector"  is  the  model 
number  of  the  detector  with  the  distance  to  the  detector  listed  in  parentheses  (the 
variations  are  due  to  slight  variations  in  the  distance  from  the  kapton  window  to  the 
detector  window  referred  to  as  the  air  gap),  "Nominal  Angle"  is  the  crystal  angle  with 
respect  to  the  electron  beam  line. 

Data  was  obtained  on  May  5^  by  rotating  the  crystal  157  5  degrees  clockwise  (G^  = 
157,5  degrees)  which  corresponds  to  Gg  =  22.5  degrees.  The  window  to  detector 
distance  (air  gap)  was  1.4  cm  of  air.  The  ORTEC  450  amplifier  was  used.  The  PXR 
spectrum  is  shown  in  Figure  9. 

Data  was  obtained  on  May  8^  at  the  same  angle  for  a  longer  time  to  improve  signal  to 
noise.  The  ORTEC  450  amplifier  was  replaced  by  the  TENNELEC  TC  244  amplifier  for 
this  and  all  the  subsequent  single  angle  experiments.  The  amplifier  was  set  for  a  12  fts 
peaking  time,  which  gives  the  best  energy  resolution.  The  delay  time  of  the  linear  gate 
was  adjusted  based  upon  observation  of  the  copper  foil  fluorescence  to  30.6  ^s.  The  air 
gap  was  l  .l  cm.  The  PXR  and  calibration  spectra  are  shown  in  Figures  10  &  1 1  The 
calibration  for  May  5^  is  similar  to  that  of  May  8*^  and  is  not  shown. 

Data  was  also  obtained  on  July  7^  at  a  Bragg  angle  of  22.5  degrees.  The  LIN  AC  was 
tuned  to  88.0  MeV  and  the  air  gap  thickness  was  1.0  cm.  The  amplifier  and  gating  setups 
were  the  same  as  those  on  May  8*.  There  appeared  to  be  a  drift  problem,  either  a  drift  in 
klystron  noise  as  the  machine  warmed  up,  or  a  detector  drift.  To  reduce  the  effects  of 
this  drift  the  experiment  was  performed  as  follows.  First  a  calibration  was  obtained  using 
the  sandwich  foil  starting  at  0  hours,  then  using  the  carbon  crystal,  a  PXR  spectrum  was 
obtained  for  two  hours.  Due  to  a  setup  error  the  crystal  was  positioned  such  that  the 
electron  beam  missed  the  tin  backing,  and  therefore  there  wasn't  a  tin  peak  for  the 
calibration  and  later  for  an  absolute  photons  per  electron  calculation  This  spectrum  was 
saved  and  a  second  calibration  was  obtained  after  three  hours.  This  was  followed  by  two 
separate  two  hour  PXR  spectra,  which  were  obtained  without  moving  the  ladder  This 
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Energy  in  KeV 


Figure  9.  PXR  spectrum  from  the  {002}  planes  of  the  carbon  graphite  crystal,  obtained 
on  May  5  at  an  angle  of  22.5  degrees.  This  spectrum  clearly  shows  6  orders.  The  low 
energy  peak  is  LINAC  induced  noise. 
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Figure  10.  PXR  spectrum  from  the  {002}  planes  of  the  carbon  graphite  crystal,  obtained 
on  May  8  at  an  angle  of  22.5  degrees.  This  spectrum  clearly  shows  7  orders.  The  pulse 
height  analyzer  discriminator  was  set  at  approximately  4  KeV  to  discriminate  out  the 
LINAC  induced  noise. 


25 


0  5  1  0  1  5  20  25  30  35  40 

Energy  11  KeV 

Figure  1 1 .  Energy  calibration  spectrum  obtained  on  May  8.  This  calibration  used  the 
spectrum  of  three  separate  foils  which  were  added  together  by  use  of  a  spreadsheet  on 
QPRO.  The  pulse  height  analyzer  discriminator  was  set  at  approximately  4  KeV  to 
discriminate  out  the  LINAC  induced  noise. 
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resulted  in  three  independent  PXR  spectra  at  a  nominal  Bragg  angle  of  22  5  degrees  The 
crystal  was  then  rotated  to  a  nominal  Bragg  angle  of  23  degrees  and  a  spectrum  was 
obtained  for  two  hours,  followed  by  a  third  calibration  after  9  hours  from  the  beginning  ot 
this  day's  data  accumulation. 

The  calibration  information  is  shown  in  Table  2.  This  table  shows  that  the  second 
calibration's  centroids  appear  to  shift  about  six  channels  to  the  left  compared  to  their 
corresponding  peaks  in  the  first  calibration,  while  the  third  calibration's  centroids  appear 
within  one  channel  of  the  corresponding  peaks  in  the  first  calibration.  This  information 
shows  that  it  is  crucial  to  have  tin  or  another  fluorescent  foil  in  place  for  calibration  while 
taking  a  PXR  spectrum.  It  also  suggests  that  data  be  taken  in  several  shorter  time  periods 
rather  then  one  longer  time  period.  With  knowledge  of  the  tin  peak's  position  for  each  run 
one  can  then  shift  the  spectra  appropriately  and  then  add  them  together  This  procedure 
will  markedly  improve  the  spectral  resolution. 

The  spectra  taken  at  22.5  degrees  were  combined  in  the  following  manner.  The  centroid 
of  the  tin  calibration  peak  in  each  spectra  was  determined.  Then  by  shifting  the  channels  of 
the  three  spectra  so  that  the  tin  centroids  were  the  same,  the  three  spectra  were  added 
channel  by  channel  on  a  QPRO  spreadsheet  [Ref  22].  Figures  12-14  show  the  combined 
PXR  spectrum  at  22.5  degrees,  the  spectrum  at  23  degrees  and  the  energy  calibration 
spectrum  respectively. 

On  July  8^  similar  data  was  obtained,  using  the  new  linear  gate  set  at  35  ps  wide  and 
delayed  from  the  machine  start  cycle  time,  Tq  by  7.7  ps.  The  air  gap  thickness  was  11 
cm.  Additionally  the  amplifier  peaking  time  was  reduced  from  12  ps  to  8  ps.  A  lower 
peaking  time  reduces  the  effect  of  double  counting  single  events  since  the  effects  of  double 
counting  are  proportional  to  the  pulse  peaking  time,  however  there  is  a  tradeoff  with 
energy  resolution.  The  energy  resolution  is  related  to  the  pulse  peaking  time.  There  are 
two  contributions.  One  is  from  parallel  noise  (detector  leakage  current,  and  all  resistors 
which  parallel  the  detector)  that  is  proportion  to  the  square  root  of  the  peaking  time.  The 
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TABLE  2  JULY  7^  CALIBRATION  DATn,  SHOWING  PEAK  ENERGY. 


CENTROID  CHANNEL  NUMBER  AND  PEAK  FWHM 


Cal. 

0  hours 

Theoretical 
Centroid 
Energy  (keV) 

Centroid 
channel  U 
(shift) 

Full  Width 
Half  Max. 
(keV) 

Titanium 

Ka 

4.509 

191 

0.640 

Yttrium 

Koc 

14932 

747 

0.620 

Kp 

16.7 

847 

0675 

Tin 

Ka 

25  196 

1307 

0,735 

Kp 

28.486 

I  1487 

1.470 

Cal 

3  hours 


Titanium 

Ka 

4.509 

■RiiH 

0.320 

-150 

Yttrium 

Ka 

14.932 

739 
_  .  .  (-7) 

0.415 

-230 

Kp 

16.7 

840 

_ iJ) _ 

0.490 

-140 

Tin 

Ka 

25,196 

0.490 

-90 

Kp 

28,486 

i 

1 

1488 

(-1) 

1.22 

-50 

Cal 

9  hours 

- j 

Titanium 

Ka 

4.509 

mmSStm 

0.350 

-70 

yttrium 

Ka 

14.932 

0.305 

1 

'-J 

o 

Kp 

16,7 

HEhI 

0370 

0 

Tin 

Ka 

25  196 

0.345 

0 

Kp 

28.486 

1488 

(^0 

0.290 

+60 

28 


Figure  12.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  July  7* 
at  an  angle  of  22.5  degrees.  This  spectrum  was  obtained  by  combining  three  separate 
shorter  time  period  PXR  spectrum.  This  spectrum  shows  six  orders  with  the  fifth  order 
and  the  tin  calibration  peak  slightly  overlapping,  with  the  tin  peak  at  a  higher  energy  than 
the  n=5  peak. 
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Figure  13.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  July  7^ 
at  an  angle  of  23.  degrees.  This  spectrum  shows  six  orders  with  the  fifth  order  and  the 
tin  calibration  peak  slightly  overlapping,  with  the  tin  peak  at  a  higher  energy  than  the  n=5 
peak. 
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Figure  14.  First  of  three  energy  calibrations  obtained  on  July  7*.  This  calibration  used 
the  sandwich  foil  containing  Titanium,  Yttrium  and  Tin.  The  and  Kp  peaks  of  Y  and 

Sn  are  both  visible,  while  the  Ti  only  shows  the  peak. 
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second  is  serial  noise  (total  input  capacitance,  transconductance  (gain  of  the  input  device 
etc.) )  that  is  inversely  proportional  to  the  square  root  of  the  peaking  time  This  results  in 
an  optimum  peaking  time  for  SiLi  detectors  between  8  and  24  |is,  [Ref  20]  The  LINAC 
was  tuned  to  84.5  MeV  Additionally,  the  PHA  was  updated  to  the  PHA-II  software 
which  has  8192  channels  compared  to  2048  for  PHA-I.  The  PHA-II  digital  stabilization 
was  utilized  to  counter  the  effects  of  drift.  An  energy  calibration  was  obtained  before  and 
after  the  PXR  spectrum.  To  compare  the  results  obtained  with  the  results  from  the  2048 
channel  analyzer  the  data  was  binned  using  a  simple  program  written  in  C  (see  Appendix 
B).  Four  channels  were  added  together  and  placed  in  one  bin  using  the  energy  value 
obtained  from  the  third  of  the  four  bins.  Figures  15  &  16  show  the  binned  PXR  spectrum 
and  the  first  binned  energy  calibration  spectnim.  The  resolution  is  probably  improved 
because  of  collecting  data  over  more  channels  and  then  binning,  and  due  to  the  use  of  the 
PHA-II  digital  stabilization  option. 

B.  SINGLE  ANGLE  PXR  ANALYSIS  (FAR  REGION) 

In  order  to  compare  the  data  with  theory,  and  from  one  experiment  to  the  next  the 
following  analysis  was  performed.  First  the  net  area  of  each  spectral  peak  was 
determined.  This  was  done  two  ways. 

The  first  way  used  the  software  supplied  with  the  PCA I  TENNELEC  PHA.  Using  the 
PHA  computer  analysis  the  peaks  were  determined  by  estimating  where  each  peak  began 
and  ended,  and  setting  these  channels  as  a  region  of  interest  The  PHA  analyzer 
determines  the  average  background  by  averaging  the  start  channel  with  the  three  previous 
channels  and  the  stop  channel  with  the  three  following  channels  and  drawing  a  straight  line 
connecting  these  two  average  values.  All  counts  below  this  noise  line  are  treated  as  noise 
and  counts  above  it  as  signal.  The  program  determines  the  gross  and  net  area.  Gross  area 
is  the  total  number  of  counts  in  the  region  (including  noise)  and  net  area  is  only  the 
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Figure  15.  PXR  spectrum  from  the  (002}  planes  of  carbon  graphite,  obtained  on  July  8* 
at  an  angle  of  22.5  degrees,  using  the  Canberra  (SiLi)  detector.  This  spectrum  shows  six 
orders  with  the  fifth  order  and  the  tin  calibration  peak  slightly  overlapping,  with  the  tin 
peak  at  a  higher  energy  than  the  n=5  peak. 
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Figure  16.  Energy  calibration  obtained  on  July  8***  This  calibration  used  the  sandwich 
foil  and  the  Canberra  (SiLi)  detector.  This  spectrum  shows  the  peak  of  Ti,  Y  and  Sn, 

and  the  Kp  peaks  of  Y  and  Sn.  This  spectrum  was  plotted  by  binning  four  channels 
together  and  plotted  at  the  central  channel  of  these  four  bins. 
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number  of  counts  above  the  noise  line  It  also  determines  the  FWHM  and  the  centroid  of 
the  peak  using  the  following  equation: 


chn 


cird 


^countS;  *chn, 
counts, 


(31) 


where  Chnj  is  the  channel  number  (or  equivalent  energy)  and  countsj  is  the  corresponding 
number  of  counts  in  that  channel. 

The  second  way  was  to  determine  the  average  noise  in  the  regions  between  the  peaks 
using  the  gross  area  of  100  channels  divided  by  100.  The  background  was  determined  by 
multiplying  the  number  of  channels  in  the  peak  by  the  average  noise  in  the  region  just 
following  the  peak.  The  net  area  is  determined  by  subtracting  this  background  from  the 
Gross  area  The  error  in  the  Net  area  is 


- %  = - ,  (32) 

since  the  error  in  the  area  is  VArea .  Tables  3-7  show  these  values  for  the  various 
experiments  Figure  17  shows  a  plot  of  the  comparison  between  the  PHA  net  area  and  the 
Extended  Background  net  area  for  the  May  8th  data.  This  plot  shows  that  except  for  the 
first  two  peaks  there  is  very  little  difference  between  the  net  areas  as  determined  by  these 
two  methods. 

The  raw  spectrum  area  needs  to  be  corrected  for  attenuation  as  the  photon  travels 
between  the  target  and  the  detector.  Also  the  relative  efficiency  of  the  detector  must  be 
taken  into  account.  The  attenuation  coefficients  were  determined  using  a  computer 
program  called  XCOM  [Ref  21]  near  the  centroid  energy  Some  of  the  coefficients  are 
listed  in  Table  8,  along  with  the  detector  efficiency  [Ref  5]  and  the  values  of  the 
correction  factors  for  Be,  Kapton  and  for  air  at  the  two  extreme  thickness  values  The 
final  column  shows  the  total  correction  factor  for  an  air  thickness  of  1 . 1  cm. 
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TABLE  3:  MOSAIC  CARBON  CRYSTAL {002}  PLANE  PEAK  DATA  TAKEN  MAY 
5th  for  90  MEV  ELECTRON  BEAM  AT  BRAGG  ANGLE  OF  0b  =  22  5  DEGREES 
FOR  PEAKS  1-6 


20  04 


4  48 


2 

3 

9.75 

14  85 

2.18 

3  32 

5 

6 

25.1 

30  3 

5  62 

6  77 

mmm 

423 

482 

202 

76 

36 

15 

5 

5 

5 

5 

5 

5 

net 

418 

477 

197 

71 

31 

10 

■■■ 

— 

13822 

— 

16570 

7261 

3286 

— 

1284 

630 

PHA 

858 

1 

588 

736 

642 

462 

320 

r  jt 

12964 

15982 

6525 

2644 

822 

310 

AN/N% 

0.93 

0.82 

1.37 

1.75 

5  1 

10 

Extended  Background 

400 

485 

455 

500 

380 

400 

net 

13422 

16525 

6806 

2786 

904 

240 

AN/N% 

0.89 

0,81 

1.29 

2  21 

4.1 

13 

TABLE  4:  MOSAIC  CARBON  CRYSTAL{002}  PLANE  PEAK  DATA  TAKEN  MAY 
gth  for  90  MEV  ELECTRON  BEAM  AT  BRAGG  ANGLE  OF  0b  =  22  5  DEGREES 


FOR  PEAKS  1-^ 

\. 

peak  # 

1 

2 

3 

4 

E(n)  keV 

4,905 

9.509 

14.239 

190301 

E(n)/E(l) _ 

1 

1.94 

2.9 

3  88 

Peak 

sross 

1213 

1301 

552 

260 

background 

8 

8 

8 

8 

net 

1205 

1293 

544 

252 

Area 

gross 

25888 

31844 

16189 

8192 

PHA 

background 

918 

869 

750 

850 

net 

24970 

30975 

15439 

7342 

AN/N% 

0.66 

0,56 

0.84 

1.3 

Extended  Background 

background 

424 

624 

592 

672 

net 

25464 

31220 

15597 

7520 

AN/N% 

0.64 

0.58 

0.83 

1.25 
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TABLES.  SAN 

DE  AS  TABLE  4  FOR  PEAKS  5-8 

iinHi 

5 

6 

7 

8 

23.828 

28  5723 

33  3524 

38  368 

E(n)/E(l) 

4.86 

5  83 

68 

7  82 

Peak 

gross 

121 

69 

33 

20 

7 

6 

5 

4 

net 

114 

63 

28 

16 

Area 

f^m 

4473 

2287 

1248 

733 

PHA 

624 

574 

352 

252 

net 

3849 

1713 

896 

481 

AN/N% 

1.85 

3.12 

4.46 

6.52 

Extended  Background 

721 

486 

435 

332 

net 

3752 

1801 

813 

401 

AN/N% 

1  92 

2.92 

5.05 

8.14 
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TABLE  6:  MOSAIC  CARBON  CRYSTAL{002}  PLANE  COMBINED  PEAK  DATA 
TAKEN  JULY  7th  FOR  88  MEV  ELECTRON  BEAM  AT  BRAGG  ANGLE  OF  = 
22,5  DEGREES  FOR  PEAKS  1-5  ANALYZED  BY  COMBINING  THE  THREE 
INDIVIDUAL  22  5  DEGREE  SPECTRU^.1 


■■1 

1 

2 

3 

4 

5 

4.85 

9  50 

1430 

19  10 

23  99 

HHniii 

1 

1.96 

2.95 

3.94 

4.95 

Peak 

mm 

4905 

- 1 

740 

287 

123 

55 

17 

17 

15 

13 

11 

net 

888 

723 

272 

110 

44 

Area 

mm 

23812 

23590 

10640 

5304 

2454 

PHA 

1315 

1269 

1517 

1434 

1213 

net 

22497 

22321 

9123 

3870 

1241 

AN/N% 

0,70 

0.71 

1.21 

2.21 

4.88 

Extended  Background 

1283 

1518 

1345 

1388 

1016 

net 

22529 

22072 

9295 

3916 

1438 

AN/N% 

0.70 

0.72 

1.18 

2.09 

4.10 
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TABLE  7.  MOSAIC  CARBON  CRYSTAL{002}  PLANE  PEAK  DATA  TAKEN  KJLY 
7th  for  88  MEV  ELECTRON  BEAM  AT  BRAGG  ANGLE  OF  0b  =  23  0  DEGREES 
FOR  PEAKS  1-5. 


■HBi 

1 

2 

3 

4 

5 

All 

9.47 

14.39 

19.32 

2420 

HRHii 

_ \ _ 

2.01 

3.05 

4.09 

5.13 

Peak 

mm 

735 

843 

349 

149 

70 

11 

10 

10 

8 

7 

net 

llA 

833 

•  339 

141 

63 

Area 

■■■ 

19766 

25832 

12319 

5766 

PHA 

1020 

— 

960 

1248 

855 

net 

18746 

24872 

11071 

4911 

22171 

AN/N% 

0.77 

0.66 

1.05 

1.66 

Extended  Background 

921 

949 

986 

795 

net 

18845 

24883 

11333 

4971 

AN/N% 

0.76 

0.66 

1.02 

1.35 

• 

^The  net  area  was  determined  by  use  of  a  peak  fitting  program  with  the  following  parameters:  Gaussian 
Peak;  25. 197  KeV  centroid,  0.548  FWHM  with  a  background  of  counts  =  20.94  -  0.527*Energy,  which  is 
consistent  with  the  balue  of  7  used  in  the  Peak  channel. 
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Not  Aroai  (count 
(Thousands) 


Figure  17.  Plot  of  the  net  area  in  each  peak  as  a  function  of  the  peak  number  as 
determined  by  the  pulse  height  analyzer  and  by  using  the  extended  background  method  for 
the  May  8th  data. 


TABLE  8.  ATTENUATION  COEFFICIENTS  FOR  THE  PEAK  ENERGIES. 
DENSITY  AND  THICKNESS  FOR  BE.  DRY  AIR  AND  KAPTON _ 


Energy 

(KeV) 

H  (cm2/g) 

Be 

H(cm2/g) 

Dry  Air 

l€  ^^^t=1.0cinl 

(1  (cm2/g) 

Kapton 

Detector 

Efficiency 

Total 
correction 
factor  for 
t  air  =1.1 

cm 

5 

wmoMm 

mSEBU 

1.0 

0.833 

10 

0.647 

rO.9941 

5.12 

10,9931  (0.9911 

3.18 

(0.9891 

1.0 

0.976 

15 

1.61 

(0.9981  (0.9971 

1.04 

(0.9961 

1.0 

0.991 

20 

0.778 

(0.9991  (0.9991 

0.542 

(0.9981 

0.97 

0,965 

25 

0.484 

(0.9991  (0.9991 

0.366 

(0.9991 

0.86 

0,857 

30 

0.179 

r0.9981 

0.354 

(1.001  (0.9991 

0.288 

(0.9991 

0.72 

0  717 

iBI 

1.82 

0.001294 

1.42 

- 

- 

thickness 

0.005  cm 

1.0<t<1.4  cm 

0.0025  cm 

- 

Using  the  net  area  from  the  PHA  analysis  and  the  coefficients  listed  above  the  corrected 
counts  is  determined  using  the  following: 

I,  =  (33) 

where  I  is  either  Net  Area  or  Net  Peak  counts  counted  by  the  (SiLi)  detector. 

Two  methods  of  analysis  were  used  to  determine  the  counts  at  the  front  face  of  the 
crystal  Iq.  The  first  and  simplest  used  equation  33  applied  peak  to  peak  where  (|ipt)total 
is  assumed  constant  over  the  entire  peak.  This  is  the  analysis  used  in  references  4  &  5 
However,  for  the  low  energy  peaks  p  varies  over  the  peak  and  the  values  obtained  for  Iq 
vary  depending  on  if  you  use  p(5  keV)  or  for  each  individual  channel.  For  this  reason  a 
second  method  was  examined.  The  second  method  involved  using  a  least  squares  fit  of 
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the  attenuation  coefficient.  The  actual  coefficients  and  the  approximations  for  the  three 
materials  are  shown  in  Figures  18-20.  The  formulas  give  |X  as  a  fimction  of  energy  and 
are  applied  to  the  values  of  energy  and  the  number  of  counts  on  a  QuattroPro  (QPRO) 
spreadsheet.  [Ref  22]  The  counts  are  then  summed  over  the  peak  and  a  new  Iq  is 
obtained.  See  Table  9  for  the  areas  and  the  ratio  of  the  area  of  the  nth  order  to  the  area  of 
the  first  order  for  the  raw  data  as  collected  by  the  detector,  the  corrected  data  using  the 
linear  attenuation  coefficient  at  the  centroid  and  the  data  corrected  using  the  linear 
attenuation  coefficient  for  each  channel.  Figure  21  plots  the  raw  area,  the  area  corrected 
for  attenuation  based  on  the  coefficient  at  the  peak  centroid  and  the  coefficient  for  each 
channel  using  the  least  square  fit  values  for  the  May  8th  data. 

It  is  not  clear  that  the  use  of  the  channel  by  channel  correction  is  appropriate  since  the 
detector  response  not  only  broadens  the  PXR  spectrum,  there  is  a  folding  of  the  detector 
response  into  the  PXR  spectral  response.  If  the  detector  resolution  is  small  compared  to 
the  PXR  resolution  the  channel  by  channel  correction  method  is  a  better  approximation,  if 
the  PXR  resolution  is  small  compared  to  the  detector  resolution  the  correction  at  the 
centroid  is  the  better  approximation,  and  if  the  two  are  comparable  which  is  the  case  for 
these  experiments  as  seen  below  the  best  method  would  be  to  deconvolve  the  PXR 
spectrum  from  the  detector  response  and  then  apply  a  channel  by  channel  correction. 
Since  this  method  is  extremely  complex  it  is  not  applied  here,  and  for  ease  of  calculation 
the  method  of  using  the  attenuation  coefficient  at  the  peak  centroid  will  be  utilized 
for  the  remainder  of  the  analysis.  In  order  to  determine  the  best  method  the  PXR 
resolution  needs  to  be  calculated. 

1.  Detector  Resolution  at  PXR  Spectral  Energies. 

To  determine  the  detector  energy  resolution,  the  calibration  peaks  were  assumed  to 
be  Gaussian  in  shape.  The  natural  line  widths  of  the  . x-ray  fluoresence  peaks  were 
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Figure  18.  Plot  of  the  log  of  the  actual  linear  attenuation  coefficient  for  Be  vs.  energy  in 
KeV,  and  the  formula  values  using  3  section  least  square  fit. 


Figure  19.  Plot  of  the  log  of  the  actual  linear  attenuation  coefficient  for  Dry  Air  vs.  energy 
in  KeV,  and  the  formula  values  using  3  section  least  square  fit. 
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Figure  20  Plot  of  the  log  of  the  actual  linear  attenuation  coefficient  for  Kapton  vs.  energy 
in  KeV,  and  the  formula  values  using  3  section  least  square  fit. 

obtained  from  Reference  23,  these  natural  line  widths  varied  from  10  eV  to  45  eV.  The 
natural  line  widths  were  also  assumed  to  be  Gaussian  distributions.  F^’orn  the  measured 
fluorescence  peaks  the  detector  resolution  or  linewidths  were  determined  by  adding  them 
in  quadrature; 

+  (34) 

where  rot,  is  the  measured  linewidth,  is  the  detector  resolution  and  F,^  is  the  natural 
fluorescent  linewidth.  The  units  are  keV.  The  uncertainty  in  F,,^  is  assumed  to  be  zero, 
therefore  the  uncertainty  in  is  the  same  as  the  uncertainty  in  F,^,  The  experimental 
values  of  Fj„  and  F^^,  by  less  than  two  tenths  of  a  percent.  These  results  are  in 
agreement  with  those  previouly  obtained  [Ref  4.]  Thus  the  value  of  the  detector 
resolution  to  be  used  will  be  that  of  the  observed  calibration  spectrum  line  widths. 
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TABLE  9  DATA  CORRECTED  FOR  ATTENUATION  THROUGH  KAPTON 
WINDOW,  AIR  AND  BE  DETECTOR  WINDOW 


I 


DiaiLvlKnvJc* 


rea  Corrected  for  attenuation  at  the  peak  centroid:  ratio  Arean/Area 


15886;  1 


30124;  1 


27000 ;  1 


15575;  0.52  7379;  0,24 


SEfi 


1246;  0.05 


25471  ;  1.13 

11 166;  0.50 

4936 ;  0,22 

29262 ;  0.85 

13786 ;  0.40 

5678;  0,16 

rea  corrected  at  each  channel  using  fitted  coefTicients;  ratio  Arean/Area 


0 

4-1 

c 

3 

0 

0 

(C 

0 


0 

z 


0 

■D 

c 

(0 

a 

3 

0 

£ 

h 


3 

Order 


4 


Raw  ^PHA 


Fitted 


Figure  21.  Plot  of  the  raw  area  at  the  detector,  the  corrected  area  using  the  linear 
attenuation  coefficient  at  the  centroid  of  the  peak  (labeled  centroid)  and  the  linear 
attenuation  coefficient  for  each  channel  based  on  a  least  square  fit  for  the  May  8th  data. 
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Table  10  lists  the  values  of  the  natural  x-ray  fluorescence  linewidths  and  the  observed 
linewidths  as  determined  by  the  pulse  height  analyzer  for  the  calibration  foils  for  the  May 
5^,  May  8^,  July  7^h  ^nd  July  spectra. 

TABLE  10.  LIST  OF  CALIBRATION  FOILS  AND  THE  RESPECTIVE  X-RAY 
ENERGY  PEAK  VALUES,  INCLUDED  IN  THE  LIST  ARE  THE  NATURAL 


LINEWIDTHS  OF  THE  FOILS,  AND  THE  OB 

SERVED  RESC 

ILUTION 

Date 

Element 

Energy 

(keV) 

I'lw 

^obs 

Cu 

8.04 

0.003 

0.698 

Y 

14.92 

0.0062 

0.628 

Cu 

8.04 

0003 

0.394 

Y 

14.92 

0.0062 

0.269 

In 

24.10 

0.0103 

0.466 

HHIH 

Ti 

4.51 

0.002 

0.325« 

Y 

14.92 

0.0062 

0.417^ 

Sn 

25.196 

0.0103 

0.493‘» 

Ti 

4.51 

0.002 

0.180 

Y 

14.92 

0.006' 

0.300 

Sn 

25.196 

0.0103 

0  248 

^The  linewidth  of  Y  wasn't  listed  in  rcferenc?  23.  The  values  was  approximated  by  using  the  values  of  Sr 
(5.6  eV)  and  Zr  (6.7  eV). 

^The  linewidth  of  In  wasn't  listed  in  reference  23,  The  value  was  approximated  by  using  the  values  of  Pd 
(10.2eV)  and  Ag  (10,0  eV). 

^The  observed  resolution  was  taken  from  the  calibration  taken  after  three  hours. 

^The  linewidth  of  Sn  wasn't  listed  in  reference  23.  The  value  was  approximated  by  using  the  values  of  Pd 
&  Ag  (10  eV).  The  next  element  listed  was  W  (43  eV)  giving  an  upper  limit  which  is  still  insignificant 
compared  the  the  observed  resolution. 
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2.  PXR  Bandwidth  Determination. 


The  PXR  bandwidth  is  determined  in  the  same  manner,  that  is 

r-  -  r'-  _r- 

*  PXR  “  *  »(»  ‘  d«i  ’ 


(35) 


where  TpxR  is  the  bandwith  of  the  PXR  spectral  line.  The  units  are  in  keV  Table  1 1  lists 
the  values  of  the  observed  resolution  for  a  given  energy,  the  calculated  PXR  bandwidth 
for  the  various  single  angle  experiments  and  the  expected  value  as  calculated  using 
equation  (27). 

Note  that  the  observed  linewidths  in  Table  10  vary  with  the  various  experiments 
These  variation  are  due  to  the  changes  made  in  the  experimental  procedure  That  is  the 
May  5th  data  is  much  wider  due  to  poorer  resolution  of  the  ORTEC  450  amplifier  This 
amplifier  was  replaced  by  the  TENNELEC  TC  244  amplifier  following  that  experiment 
The  July  8th  data  used  an  upgraded  version  of  the  PHA  software  which  has  digital 
stabilization  features  which  compensate  for  drift  and  therefore  reduces  the  observed 
linewidth.  Notice  that  in  spite  of  these  changes  in  detector  resolution,  the  PXR  linewidths 
do  not  vary  significantly,  since  the  changes  in  the  detector  resolution  are  accounted  for  by 
equation  (35). 

TABLE  11  LIST  OF  OBSERVED  RESOLUTION  AND  CALCULATED  PXR 
BANDWIDTH  FOR  THE  VARIOUS  SINGLE  ANGLE  EXPERIMENTS 
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TABLE  11.  cont. 


Date 

Energy  (keV) 

__  Tobs  (keV) 

WSIMSSSM 

Theory 

7  July  22.5 
degrees 

4  8 

047^ 

0  34 

0  17 

9  8 

0.54^ 

0.39 

0  35 

14  3 

063'* 

0  38 

0  52 

19  1 

0  54-‘ 

0.30 

0  70 

24.0 

1.02’ 

090 

0  87 

5 

0.47 

0.34 

0  17 

15 

0.61 

042 

0.50 

25 

0.82 

0.65 

0.83 

5 

0.37 

0.17 

15 

0.52 

0.43 

0.52 

25 

0.26^ 

020 

087 

Once  the  resolution  and  the  net  area  is  determined  at  the  front  face  of  the  crystal, 
these  values  can  be  compared  to  the  theoretical  values.  Since,  for  most  of  these 
experiments,  there  was  no  reference  tin  peak  from  which  to  compute  the  LINAC  current, 
the  analysis  is  accomplished  by  comparing  the  experimental  ratios  of  the  intensity  (area)  of 
a  given  order  to  the  intensity  of  the  first  order  with  that  calculated  from  theory  The  single 
angle  experiments  were  all  obtained  in  the  far  region,  therefore  the  theoretical  ratios  are 
determined  by  integrating  equation  (26)  over  the  field  of  view,  and  the  bandwidth  is  based 
on  Ad^  and  not  the  FWHM  of  J2  The  experimental  ratios  are  listed  in  Table  9  above.  As 
described  in  Chapter  II  ,  section  C,  the  calculations  must  be  performed  over  the  limited 
field  of  view. 

Figure  2  in  Chapter  II  section  C  illustrates  how  the  aperture  in  the  far  case  is  less 
than  the  FWHM  of  J2  The  theoretical  ratios  are  calculated  by  integrating  equation  (26) 
over  the  field  of  view  for  each  order  then  taking  the  ratio  of  the  number  for  a  given  order 
to  the  first  order  Since  the  effects  of  mosaic  spread  and  multiple  scattering  aren't  well 


*The  observed  resolution  was  taken  from  the  calibration  taken  after  three  hours. 
^Estimated  due  to  overlap  of  fifth  order  peak  with  the  tin  peak  in  the  PXR  spectrum. 
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understood,  the  theoretical  values  with  and  without  these  effects  are  included  for 
comparison.  These  ratios  are  shown  in  Table  12  along  with  the  experimental  ratios  The 
experimental  ratios  are  plotted  as  a  function  of  order  in  Figure  22.  The  experimental 
values  are  compared  with  the  theory  by  dividing  the  experimental  value  by  the  theoretical 
value  and  is  also  shown  in  Table  12 


TABLE  12  COMPARISION  OF  EXPERIMENTAL  AREA  RATIOS,  I(n)/I(l).  WITH 
THEORY  WITH  AND  WITHOUT  THE  EFFECTS  OF  MULTIPLE  SCATTERING 
AND  MOSAIC  SPREAD 


Experimental  Data 


15900 


16400 


6600 


2700 


800 


30200 


31700 


15600 


7400 


3900 


27000 


22900 


9200 


3900 


1200 


22500 


25500 


11200 


4900 


2200 


34600 


29300 


13800 


5700 


1300 


Ratio 


Theory  without 
scattering  or  mosaic 
spread 


Ratio 


0.27 


0.05 


0.01 


0.005 


0.02 


0.004 


0.001 


0.0003 


1 


0.27 


0.05 


0.01 


0.005 


Exp 

theory 


1 


3  8 


8.4 


17 


10 


Theory  with 
scattering  and 
mosaic  spread 


Ratio 


0.08 


0.01 


0.003 


0.001 


Figure  22.  Plot  of  Area  ratio  of  nth  order  peak  to  first  order  peak  vs  order  for  the  four 
experiments  conducted  near  the  Bragg  angle  of  22.S  degreees. 


Notice  that  the  experimental  ratios  are  consistantly  larger  than  theory,  and  the  inclusion  of 
mosaic  spread  and  multiple  scattering  makes  the  discrepancy  even  larger,  thus  showing 
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that  equation  (21)  may  be  preferred  over  equation  (20)  This  indicates  that  the  theory  is 
inadequate  to  explain  the  observed  results  It  also  indicates  that  it  is  important  to  look  at 
the  convolution  of  the  electron  beam  divergence  with  the  multiple  scattering  effect  for  an 
even  better  estimation  of  the  effect  of  multiple  scattering. 

3.  Calculations  of  Absolute  PXR  Yield  (Photons/Electron) 

In  addition  to  comparing  the  intensity  ratios,  for  the  7  and  8  July  experiments  the 
number  of  photons  per  electron  can  be  calculated  The  Linac  current  can  be  determined 
by  using  the  yield  of  the  Tin  fluorescent  peak  using  the  following  equation 


^  Area 

V 


'  LWAC 


^Tin 


(36) 


where  IlinaC  **  calculated  beam  current,  Areajjj,  is  the  Area  of  the  tin  peak  at  the 
detector  (from  the  pulse  height  analyzer),  time  is  the  collection  time  of  the  PXR  spectrum, 
AW  is  the  atomic  weight,  a  is  the  electron  interaction  cross  section  (1.0  e-26  m2)[Ref 
24],  is  the  solid  angle  subtended  by  the  detector  (2xl0-'*  sr  ),  p  is  the  density  of  the 
target  (2.260  g/cm^),  is  Avogadro's  number,  fj^^  js  the  de-excitation  transition 
probability  (0.712)[Ref  25],  t  is  the  thickness  of  the  target  (2.98xl0'5m  for  22  5  degrees), 
e  is  the  relative  detector  efficiency  at  the  tin  peak  energy  (0.9  for  these  experiments)  and  a 
is  the  total  attenuation  factor  of  photons  (through  the  effective  thickness  of  tin,  the 
effective  thickness  of  the  mosaic  crystal,  the  kapton  window,  the  air  gap  and  the  beryllium 
window  of  the  detector  (0.763  for  22  5  degrees.)[Ref  21]). 

The  electron  interaction  cross  sections  [Ref  24]  were  obtained  for  85  MeV 
electrons  The  transition  probabilities  [Ref  25]  include  both  the  K-Ll  1  and  the  K-Ll  1 1 
de-excitation  probabilities.  The  attenuation  coefficie  •  ’  were  obtained  using  XCOM  at  the 
known  de-excitation  energy  of  Tin.  The  details  of  these  calculations  may  be  found  in 
reference  17. 


Using  the  LINAC  current  in  electrons  per  second,  the  intensity  of  the  PXR  in  counts 
as  listed  in  Table  9  and  the  data  collection  time,  the  number  of  photons  can  be  calculated 
as; 

fPXR  Area, 
j _  Time 

LINAC  Current 

These  values  are  shown  in  Table  13  along  with  the  calculated  currents. 


TABLE  13  COMPARISON  OF  THE  EXPERIMENTAL  NUMBER  OF  PHOTONS 
PER  ELECTRON  TO  THEORY  WITH  AND  WITHOUT  THE  EFFECTS  OF  MOSAIC 
SPREAD  AND  MULTIPLE  SCATTERING. 


N 


EXP 

DATA 

1 

THEORY 

WITH 

SCATT 

EXP 

THEORY 

THEORY 

WITHOUT 

SCATT 

EXP 

THEORY 

1.57e-6 

3.64e-8 

43 

1.39e-7 

11 

1.48e-6 

3.14e-9 

471 

3  69e-8 

40 

5.92e-7 

4.24e-10 

1400 

7.14e-9 

83 

2.62e-7 

1.07e-10 

2450 

1.98e-9 

132 

5.6e-8 

3.56e-ll 

2120 

6.85e-10 

110 

l.60e-6 

3.03e-8 

53 

6.24e-8 

26 

1.81e-6 

3.38e-9 

536 

1.27e-8 

143 

7.94e-7 

5.03e-I0 

1580 

2.30e-9 

345 

3.48e-7 

1.29e-10 

2700 

6.24e-10 

558 

1.56e.7 

4.33e-ll 

3600 

2.13e-10 

732 

1.58e-6 

3.63e-8 

44 

1.39e-7 

11 

1.34e-6 

3,13e-9 

428 

369e-8 

36 

6.30e-7 

4.23e-10 

1490 

7.12e-9 

89 

2.60e-7 

I.06e-10 

2450 

1.98e-9 

131 

5.94e-8 

3.55e-ll 

1670 

6.83e-10 

87 
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The  number  of  photons  per  electron  are  plotted  vs  order  for  these  three  dates  in 
Figure  23.  Note  that  the  data  is  consistently  larger  than  the  theory,  as  seen  with  the  ratios 
Figure  24  is  a  plot  of  the  number  of  photons  per  electron  for  experimental  data,  and  the 
theory  with  and  without  the  effects  of  multiple  scattering  and  mosaic  spread.  With  the 
actual  number  of  photons  per  electron  there  is  more  information  than  with  the  ratios  alone. 
First  it  is  clear  that  the  data  is  much  larger  than  expected.  Second,  the  difference  between 
the  data  and  the  theory  gets  larger  with  order,  which  explains  why  the  error  in  the  ratios 
increased  with  order.  Because  of  overlap,  a  peak  fitting  routine  was  used  to  obtain  the 
area  under  the  fifth  order  PXR  and  tin  peaks.  [Ref  26] 
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Figure  23.  Plot  of  the  d)solute  number  of  photons  per  electron  vs  order  for  the  three 
experiments  near  the  Bragg  angle  of  22.S  degrees  with  a  tin  peak  used  to  determine  the 
LINAC  current. 
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Figure  24.  Plot  of  number  of  photons  per  electron  vs  angle  for  the  7  July  22.5  degree  data 
along  with  the  theory  both  with  and  without  the  effects  of  multiple  scattering  and  mosaic 
spread. 
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C  PXR  DATA  COLLECTION  (DISPLACEMENT  FROM  THE  BRAGG  ANGLE) 

Two  sets  of  "rocking  data"  wc’e  obtained.  For  these  experiments  rocking  data  refers  to 
a  fixed  detector  angle  of  45  degrees  with  respect  to  the  beam  line,  and  the  crystal  angle  is 
rocked  about  the  Bragg  angle  of  22  5  degrees.  The  first  was  obtained  on  June  5,  1992. 
Data  was  obtained  using  the  ORTEC  SiLi  detector  and  the  TENNELEC  amplifier  This  is 
the  far  region.  The  linear  gate  was  set  at  5  |is  wide  with  a  delay  time  of  30.6  ps  from  Tq 
(the  LINAC  start  sequence  time).  The  air  gap  thickness  was  1.1  cm.  Brief  spectra  were 
obtained  starting  at  22.5  degrees  for  200  seconds.  Then  the  edges  were  determined  by 
obtaining  spectra  at  26.5  and  18.5  degrees,  both  of  which  showed  no  PXR  peaks.  To 
obtain  an  estimate  of  the  angle  giving  the  largest  count  rate  three  more  angles  were  set  for 
200  seconds  each.  The  results  are  shown  in  Table  14,  where  "counts"  is  the  net  number  of 
counts  in  the  peak  channel  of  the  n=l  peak.  The  count  rate  appeared  to  increase  upon 
rotation  from  22.5  degrees  to  lower  angles  with  the  maximum  near  22  1  degrees  Angles 
were  selected  assuming  22,1  degrees  to  be  the  angle  which  produced  the  maximum  count 
rate.  (This  doesn't  agree  with  theory  where  22.5  degrees  should  be  the  maximum,  but  due 
to  uncertainties  in  the  crystal  orientation  this  was  assumed  within  error). 


TABLE  14  NUMBER  OF  COUNTS  VS.  ANGLE  FOR  200  SECOND  SPECTRUM 


Angle 

counts 

22.5 

13 

26.5 

0 

18.5 

0 

22,3 

15 

22.1 

17 

21.9 

10 
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The  crystal  angle  was  varied  starting  at  19.1  degrees,  and  increased  by  10  degree. 
When  the  angle  was  within  1  degree  of  the  assumed  maximum  the  angle  was  increased  by 
0.5  degrees.  The  spectra  were  obtained  for  30  minutes  each.  As  seen  in  Table  15  the 
number  of  counts  in  the  peak  channel  in  the  n=l  peak  continued  to  increase  beyond  22. 1 
degrees.  The  angle  was  increased  in  0.5  degree  steps  until  the  number  of  counts  in  the 
peak  channel  were  comparable  to  the  number  at  21 . 1  degrees,  after  which  1 .0  degree  steps 
were  taken  until  very  little  n=l  peak  showed.  These  uncorrected  PXR  spectra  are  shown 
in  Figures  25  &  26.  Full  scale  spectra  can  be  found  in  Appendix  A.  This  data  shows  that 
the  maximum  count  rate  is  obtained  near  22.5  degrees  as  expected. 

On  July  22.  1992.  rocking  data  was  obtained  using  the  Canberra  (SiLi)  detector  and  the 
TENNELEC  amplifier.  This  is  the  near  region.  The  LINAC  was  tuned  to  95.9  MeV. 
The  air  gap  was  1 . 1  cm.  Due  to  a  set  up  error  the  data  collection  time  was  not  fixed  and 
therefore  varied  at  different  angles  but  was  approximately  45  minutes.  A  longer  time  was 
used  than  on  June  5^  since  the  Tin  peak  wasn't  as  well  defined.  The  angles  where  taken 
to  be  the  same  as  those  used  on  June  5^  for  comparison  purposes  with  the  exception  that 
no  spectrum  was  collected  at  19.1  degrees.  These  uncorrected  PXR  spectra  are  shown  in 
Figures  27  &  28  Full  scale  spectra  can  be  found  in  Appendix  A. 

D.  PXR  ANALYSIS  (DISPLACEMENT  FROM  THE  BRAGG  ANGLE) 

The  analysis  of  the  data  involved  three  issues.  The  first  is  the  value  of  the  PXR  energy 
versus  angle,  the  second  is  the  intensities  versus  angle  and  the  last  is  the  number  of 
photons  per  electron  versus  angle. 

To  determine  the  PXR  energy  versus  angle  the  original  spectra  were  analyzed  using  the 
pulse  height  analyzer  to  determine  the  centroid  of  each  peak,  including  the  centroid  of  the 
tin  peak.  The  value  of  the  first  order  centroid  was  plotted  vs.  angle,  along  with  the 
theoretical  value.  There  appeared  to  be  a  discrepancy  between  the  theoiy  and  the  data. 
Next  the  tin  peak  centroid  was  plotted  vs  angle.  Note  that  theoretically  the  tin  peak 
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Figure  25.  PXR  spectra  from  the  {002}  planes  of  carbon  graphite  for  the  "far  region", 
obtained  on  June  5th  at  angles  19.1-22.6  degrees. 
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position  in  the  spectrum  should  not  depend  on  angle.  The  fact  that  the  tin  peak's  position 
was  not  constant  with  angle  showed  that  there  was  some  shifting  in  energy  with  time 
Taking  the  difference  between  the  experimental  tin  peak  centroid  and  the  known  tin  de¬ 
excitation  energy  (25  1955KeV)  and  subtracting  this  difference  from  the  PXR  centroid, 
this  corrected  PXR  centroid  was  compared  with  theory  and  showed  much  better 
agreement.  These  plots  are  shown  in  Figures  29  &  30  for  the  June  5^  and  July  22"*^ 
experiments  respectively.  The  scale  on  the  left  refers  to  the  PXR  energies.  The  original 
PXR  centroids  are  plotted  as  "pxr  centroid",  the  PXR  centroid  minus  the  difference  in  the 
true  tin  and  known  tin  is  plotted  as  "corrected  cent"  The  experimental  tin  peak  is  shown 
using  the  scale  on  the  right,  with  the  known  tin  energy  plotted  as  a  dotted  line.  Note  that 
the  same  increments  were  used  on  the  two  scales  for  comparison  purposes 

This  also  explains  the  difference  between  the  expected  values  of  E(n)/E(l)  equal  to  n 
and  the  actual  values  as  listed  in  Tables  2-6.  This  is  further  illustrated  in  Table  16.  The 
(*)  indicates  corrected  values  based  on  the  shift  in  the  Tin  peak  centroid  energy  from  the 
known  de-excitation  energy.  This  shift  is  listed  as  "Tin  Shift", 


TABLE  16  COMPARISON  OF  RAW  ENERGY  RATIOS  TO  THE  RATIOS  OF  THE 
SHIFTED  ENERGIES  FOR  THE  5  JUNE  ROCKING  DATA 


ANGLE 

E(l) 

E(2) 

E(l)/E(2) 

TIN 

SHIFT 

E(l)* 

E(2)* 

E(l)*/E(2)* 

21,1 

4.73 

9.23 

1.95 

0.27 

446 

8.96 

2.01 

21.6 

4.78 

9.34 

1  95 

0.22 

4.56 

9  12 

2.00 

22,1 

4.84 

9.52 

1.97 

0.26 

4  58 

9.26 

2.02 

22.6 

4.94 

9.72 

1.97 

0.23 

4.71 

949 

2.01 

23.1 

4.98 

9.81 

1.97 

0.13 

4.85 

9.68 

200 

23.6 

5.02 

9.95 

1.98 

0.17 

4  85 

9.78 

2.02 

24,1 

5.21 

10.17 

1.95 

0.18 

5.03 

999 

1  99 

246 

5.23 

10.40 

1  99 

0  18 

5.05 

10.22 

2.02 
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Figure  30.  Plot  of  PXR  first  order  energy  original  data,  corrected  data  and  theory  vs 
angle,  tin  centroid  and  theory  for  the  July  22*^^  experiment  (near  region). 
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These  values  were  computed  by  finding  the  energy  centroid  based  on  equations  (23)  and 
(24)  within  the  limits  of  the  field  of  view.  That  is: 


where  the  integrals  are  computed  over  the  field  of  view  Ad^.  Figures  31-33  illustrate  the 
effect  of  changing  angles  on  the  field  of  view  Figure  31  illustrates  A0,  for  the  far  case 
near  the  Bragg  angle,  while  Figure  32  is  approximately  one  degree  displaced  from  the 
Bragg  angle  for  the  far  case.  Notice  that  the  effect  of  changing  the  crystal  angle  is  to  shift 
the  field  of  view  across  the  function  J2  The  dector  sees  a  slice  of  J2  and  that  slice  shifts 
across  J2  as  the  crystal  angle  changes  from  the  Bragg  angle.  Figure  33  illustrates  the  field 
of  view  for  the  near  case.  Notice  that  the  field  of  view  is  larger  than  the  FWHM  of  J2  for 
this  case. 

Comparison  of  Figures  29  &  30  show  that  for  the  far  region  the  detail  of  the  function  J2 
is  not  important  in  the  energy  centroid  calculations,  however  for  the  near  region  the  detail 
of  J2  has  a  definite  effect  on  the  calculations.  The  dip  in  Figure  30  near  23  degrees  shows 
the  effect  of  the  dip  of  J2  There  is  a  larger  discrepancy  between  the  theory  and  the 
corrected  data  for  the  July  22nd  experiment.  This  may  be  due  to  the  small  tin  peak  area 
and  therefore  a  larger  error  in  determining  the  tin  peak  centroid.  Also,  since  the  details  of 
J2  have  a  larger  effect,  an  error  in  the  angle  determination  has  a  larger  effect,  and  shifting 
the  data  by  1/4  degree  would  show  better  agreement  with  the  corrected  data  and  theory. 
Figures  34  and  35  shows  a  plot  of  the  corrected  values  of  the  energies  for  the  first  three 
orders  versus  angle,  along  with  the  theoretical  values  for  the  June  5*^  and  July  22*’*^ 
experiments  respectively.  Note  that  on  this  scale  the  difference  between  theory  and 
experiment  is  very  small,  although  there  is  a  larger  difference  at  the  smaller  angles  for  the 
higher  orders. 
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Figure  33  Illustration  of  J2  and  the  field  of  view  for  the  near  region  near  the  Bragg  angle 
of  22. 5  degrees. 


70 


Angle  in  degrees 


-*-Peaki  -♦-Peak2  -t-PeakS 
Theory  pkl  Theory  pk2  "B-  Theory  pk3 


Figure  34.  Plot  of  PXR  corrected  data  and  theory  vs.  angle  for  the  first  three  orders  for 
the  June  5^  experiment  (far  region). 
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To  compare  the  intensities  at  each  angle  the  spectra  obtained  on  June  5^  and  July  22"^^ 
at  various  angles  needed  to  be  normalized  .  Since  the  LINAC  current  was  adjusted  to 
maintain  a  machine  pulse  to  PHA  count  rate  ratio  of  about  6:1,  and  since  angles  closer  to 
the  Bragg  angle  of  22.5  degrees  produce  a  larger  number  of  x-rays,  the  current  at  these 
angles  was  lower. 

For  both  dates,  one  spectrum  was  used  to  "strip"  all  the  other  spectra  by  shifting  and 
scaling  the  spectrum  based  on  the  tin  peaks.  For  the  June  5^  data  the  19  1  degree 
spectrum  was  used  and  for  the  July  22'^**  data  the  20. 1  spectrum  was  used.  The  scaling 
factor  was  determined  by  equation  33. 


Scale  Factor  = 


Area  Tin-  Area  Ref  Tin 
Area  Tin 


(37) 


where  Area  Tin  is  the  actual  tin  area  of  the  spectrum  being  scaled  and  Area  Ref  Tin  is  the 
corresponding  area  for  the  equivalent  channels  of  the  19,1  or  20. 1  degree  spectrum  tin 
peak.  To  account  for  the  shift  in  tin  energy,  the  spectrum  was  shifted  and  the  reference 
was  subtracted  until  the  final  spectrum  showed  an  approximately  zero  net  area. 

After  the  spectrum  were  scaled  and  shifted,  the  June  5^  data  showed  an  average 
background  of  zero.  For  the  July  22"<^  data  there  remained  small  amounts  of  LINAC 
induced  noise,  especially  at  the  lower  energies.  Using  QPRO  spreadsheets  and  the  XCOM 
program  the  counts  were  corrected  channel  by  channel  for  attenuation  through  the  Kapton 
window,  the  air  gap  and  the  Be  window  of  the  detector.  The  intensity  (area)  at  the  front 
face  of  the  crystal  was  determined  by  summing  over  the  channels  of  each  peak.  The  areas 
are  listed  in  Tables  16  and  17  for  the  far  region  (June  5th  experiment)  and  the  near  region 
(July  22nd  experiment)  respectively.  The  scaled,  shifted  and  corrected  spectra  for  both 
the  June  5*^  and  July  22**^^  experiments  can  be  found  in  Appendix  A. 

Using  the  scaled  spectra,  the  intensity  can  be  plotted  as  a  function  of  angle  for  the 
various  orders.  Figures  36  and  37  show  the  scaled  area  intensities  vs.  angle  for  the  June 
5^  and  July  22"*^  experiments  respectively  Note  that  the  second  order  peak  is  larger  than 
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TABLE  16.  SCALED  AREA  INSTENSITIES  FOR  THE  JUNE  5TH  EXPERIMENT 


(FAR  REGION) 


TABLE  17  SCALED  AREA  INTENSITIES  FOR  THE  JULY  22ND  EXPERIMENT 


(NEAR  REG] 

[ON) 

ANGLE 

n=l 

2 

3 

21.1 

3850 

1860 

430 

21  6 

8680 

5970 

2120 

22.1 

13950 

11225 

5280 

22.6 

18250 

14910 

6660 

23.1 

12475 

10630 

4580 

23.6 

5040 

3130 

1200 

24.1 

2635 

1215 

24.6 

2160 

700 

25.6 


1670 


Inf  08 


Figure  36.  Plot  of  scaled  area  intensity  for  all  orders  versus  angle  for  the  June  5^ 
experiment. 
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Figure  37.  Plot  of  scaled  area  intensities  for  all  orders  versus  angle  for  the  July  22”^^ 
expen  rents. 
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the  first  order  peak  of  the  June  5^  data  at  both  22  5  and  23  degrees,  but  it  is  smaller  than 
the  first  order  peak  for  all  the  angles  of  the  July  22"*^  experiments  Also  note  that  the  plot 
as  a  function  of  angle  for  a  given  order  is  broader  for  the  July  22'^*^  experiment  than  for  the 
June  experiment  i  e  ,  the  PXR  intensity  falls  off  faster  as  a  function  of  angle  for  the  far 
region  than  for  the  near  region 

The  next  comparison  of  the  data  with  theory  involved  analysis  using  ratios  For  the  first 
order  only,  the  intensities  were  divided  by  the  first  order  intensity  at  22  6  degrees  (the 
angle  with  maximum  intensity)  Although  the  theoretical  maximum  is  at  22  5  degrees,  to 
compare  with  the  experimental  data,  the  number  of  photons  per  electron  for  each  angle 
was  divided  by  the  theoretical  value  at  23  degrees  Figures  38  &  39  show  these 
comparisons  for  the  June  5^  and  July  22"^*  experiments  respectively  Note  that  the  June 
5^  data  appears  to  be  shifted  by  approximately  0  3  degrees,  as  shown  by  the  dotted  line 
<\lso  Note  that  the  July  22nd  data  appears  to  be  shifted  in  the  other  direction  (compared  to 
the  June  5th  data)  by  about  0. 1  degrees. 

These  ratios  compared  one  angle's  first  order  intensity  to  that  of  the  first  order  intensity 
at  22  6  degrees  A  second  comparison  involves  the  ratio  of  any  given  order  intensity  to 
the  first  order  intensity  for  a  single  angle  Tables  18  and  19  show  these  theoretical  ratios 
with  the  experimental  data  in  brackets  for  the  far  and  near  case  respectively  along  with  the 
experimental  ratios  in  brackets  Figures  40  and  41  plot  the  experimental  and  theoretical 
values  for  I^2)/I(l)  and  for  I(3)/I(l)  for  both  the  far  and  near  region  respectively  Note 
that  there  is  little  variation  of  the  theoretical  ratios  with  angle  and  that  there  is  little 
difference  between  the  theoretical  values  for  the  far  and  near  case  The  fact  that  the 
experimental  values  vary  with  angle  with  a  maximum  at  9b  shows  that  searching  for  the 
maximum  in  this  ratio  can  be  used  as  a  means  to  position  the  crystal  at  65  if  the  LIN  AC 
current  isn't  monitored 

Lastly,  the  number  of  photons  per  electron  for  the  data  at  each  angle  can  be  compared 
to  the  theory  The  method  of  analysis  is  the  same  as  that  used  in  the  single  angle 
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Figure  38  Plot  of  scaled  area  intensity  divided  by  the  area  of  22.6  degrees  for  the  first 
order  versus  angle,  the  theoretical  ratio  versus  angle  and  the  shifted  theoretical  ratio  for 
the  June  5^  experiment  (far  region). 
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Figure  39  Plot  of  scaled  area  intensity  divided  by  the  area  of  22.6  degrees  for  the  first 
order  versus  angle  and  the  theoretical  ratio  versus  angle  for  the  July  22"^  experiments 
(near  region). 
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TABLE  18  RATIO  OF  INTENSITY  TO  THE  INTENSITY  OF  THE  FIRST  ORDER 
FOR  VARIOUS  ANGLES  FOR  THE  JUNE  5th  EXPERIMENT 


ANGLE 


I(2)/I(l);  [DATA]  |  I(3)/I(l);  (DATAl  I(4)/I(l);  (DATA 


II 


0034 


0.032 


0009 


0  009 


0  009 


0  009 


0008 


0  008 


TABLE  19  RATIO  OF  INTENSITY  TO  THE  INTENSITY  OF  THE  FIRST  ORDER 
FOR  VARIOUS  ANGLES  FOR  THE  JULY  22"^  EXPERIMENT 


ANGLE 


I(2)/I(l);  (DATA]  !  I{3)/I(l);  IDATA)  I(4)/I(l);  (DATA 


I 


0.034 


0  009 


0.009 


25.6 


0.182 


0.031 


0008 


I 


Figure  40.  Plot  of  the  experimental  and  theoretical  ratio  of  the  nth  order  intensities  to  the 
first  order  (I(n)/I(l))  for  the  n=2  and  n=3  peaks  for  the  June  5th  experiment  (far  region). 


n= 2  theory  A  n  =  2  data  .  n=3  theory  *  n=3data 

Figure  41 .  Plot  of  the  experimental  and  theoretical  ratio  of  the  nth  order  intensities  to  the 
first  order  (I(nyi(l))  for  the  n=2  and  n=3  peaks  for  the  July  22nd  experiment  (near 
region). 


experiments.  The  current  was  determined  by  using  the  area  of  the  tin  peak  in  the  19  1  and 
21.1  degrees  spectra  for  the  June  5*^  and  July  22”*^  experiments  respectively  These  angle 
were  used  since  the  scaling  is  based  on  them,  and  therefore  the  current  will  be  normalized 
to  this  current.  Table  20  and  21  show  the  experimental  values,  the  theoretical  values  and 
the  ratio  of  experiment  to  theory  for  the  first  order  peaks  at  the  various  angles  for  the  June 
5^  and  July  22*^*^  experiments  respectively.  The  theoretical  values  are  calculated  using  the 
theory  presented  in  Chapter  I.  The  large  variation  between  the  far  region  and  the  near 
region  are  due  to  the  larger  solid  angle  subtended  by  the  detector  in  the  near  region,  A 
larger  solid  angle  corresponds  to  a  larger  field  of  view.  These  values  were  computed 
using  computer  integration  of  equation  (23)  An  example  of  the  computer  variables  is 
presented  in  Appendix  B.  Again,  the  theoretical  values  are  much  lower  than  the  data  in 
both  the  far  and  the  near  case.  All  of  the  analysis  consistently  shows  that  the  data  is 
clearly  much  la^'ger  than  expected. 

TABLE  20.  NUMBER  OF  PHOTONS  PER  ELECTRON  FOR  THE  DATA,  THEORY 
WITHOUT  MOSAIC  SPREAD  OR  MULTIPLE  SCATTERING  AND  THEORY  WITH 
MOSAIC  SPREAD  AND  MULTIPLE  SCATTERING,  FOR  THE  FIRST  ORDER 


:AK  versus  ANGLE 

FOR  THE  JUNE  5*  EXPERII^ 

ENT  (1  =  3  39E-19  AMPS) 

ANGLE 

DATA 

THEORY 

WITHOUT 

SCATTERING 

DATA 

THEORY 

THEORY 

WITH 

SCATTERING 

DATA 

THEORY 

21.1 

8.77e-7 

8.85e-9 

10 

5  77e-9 

15 

21.6 

2.16e-7 

2.13e-8 

10 

1.30e-8 

17 

22.1 

4,64e-7 

799e-8 

508 

3.43e-8 

14 

22.6 

l,95e-6 

1.34e-7 

15 

3.77e.8 

52 

23.1 

1.65e-6 

4.72e-8 

35 

2.55e-8 

65 

23.6 

3.82e-7 

1.60e-8 

24 

1.03e-8 

37 

24.1 

1.90e-7 

7.82e-9 

24 

5.29e-9 

36 

24.6 

1.40e-7 

4.64e-9 

30 

3  20e-9 

44 

256 

5  32e-8 

2.18e-9 

24 

1  53e-9 

35 
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TABLE  21  NUMBER  OF  PHOTONS  PER  ELECTRON  FOR  THE  DATA,  THEORY 
WITHOUT  MOSAIC  SPREAD  OR  MULTIPLE  SCATTERING  AND  THEORY  WITH 
MOSAIC  SPREAD  AND  MULTIPLE  SCATTERING,  FOR  THE  FIRST  ORDER 
PEAK  VERSUS  ANGLE  FOR  THE  JULY  22"^  EXPERIMENT  (I  =  4,17X10*15 
AMPS). _ 


ANGLE 

DATA 

THEORY 

WITHOUT 

SCATTERING 

DATA 

THEORY 

THEORY 

WITH 

SCATTERING 

DATA 

THEORY 

20,1 

l,39e-6 

3  67e-8 

38 

2.44e-8 

57 

21.1 

5.26e-6 

l.32e-7 

40 

4.17e-8 

126 

21.6 

1.19e-5 

3.75e-7 

32 

8.36e-8 

142 

22.1 

1.91e-5 

8.12e-7 

24 

2.01e-7 

95 

22.6 

2.59e-5 

9.38e-7 

28 

4.33e-7 

60 

23,1 

1.63e-5 

7.31e-7 

22 

3  44e-7 

47 

23.6 

6.89e-6 

2.75e-7 

25 

1.66e-7 

42 

24.1 

3.61e-6 

1.15e-7 

31 

7.63e-8 

47 

24.6 

2.96e-6 

6.27e-8 

47 

4  30e-8 

69 

25.6 

2.28e-6 

2.76e-8 

83 

1  94e-8 

118 

Lastly,  since  the  experiment  in  the  near  region  was  done  as  "rocking  data"  there  wasn't 
a  bandwidth  resolution  calculation  computed  for  the  near  region  in  section  IV  B  1  Table 
22.  lists  the  \  ~lues  of  the  natural  linewidths  of  the  calibration  foils  and  the  observed 
resolution  for  the  July  22nd  experiment.  Table  23.  lists  the  values  of  the  observed 
resolution  for  a  given  energy,  the  calculated  PXR  bandwidth  and  the  expected  value  as 
calculated  for  the  near  region. 
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TABLE  22.  LIST  OF  CALIBRATION  FOILS  AND  THE  RESPECTIVE  X-RAY 
ENERGY  PEAK  VALUES.  INCLUDED  IN  THE  LIST  ARE  THE  NATURAL 


LINEWIDTHS  OF  THE  FOILS,  AND  TH 

E  OBSERVED  RESO 

LUTION 

Element 

Energy  (KeV) 

HERS&HI 

■IBSSH 

Ti 

4.51 

0.002 

0.150 

Y 

14.92 

0.006 

0.470 

Sn 

25.196 

0.010 

.8 

TABLE  23.  LIST  OF  OBSERVED  RESOLUTION  AND  CALCULATED  PXR 
BANDWIDTH  FOR  THE  NEAR  REGION  FROM  THE  22  1  DEGREE  PXR 


SPECTRUM  OBTAl 

"NED  JULY  22ND  (Nl 

EAR  REGION) 

(KeV) 

Tpv.  (KeV) 

Theory 

5  02 

0,450 

0.420 

0.332 

9.94 

0.500 

0.384 

0.587 

14.89 

0.205 

0.835 

*Unable  to  determine  due  to  poor  statistics  in  peak. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

Enhanced  higher  order  PXR  spectra  from  the  {002}  planes  of  carbon  have  been 
observed.  The  ratio  of  the  higher  order  peaks  to  the  first  order  peaks  (I{n)/I(l))  were 
examined  These  ratios  were  found  to  be  much  larger  than  expected,  and  the  difference 
increases  with  increasing  order  The  absolute  intensity  was  determined  based  on  the 
LINAC  current  calculations  from  the  fluorescent  x-ray  yield  of  a  0  0275  mm  thick  tin  foil 
placed  behind  the  graphite  crystal  The  absolute  intensity  was  found  to  exceed  the 
theoretical  values  by  a  minimum  factor  of  approximately  10  when  the  effect  of 
multiple  scattering  of  the  electron  beam  and  the  effect  of  the  mosaic  spread  of  the 
crystal  are  not  considered.  This  factor  increased  to  more  than  3000  when  these 
effects  are  included. 

This  is  the  first  time  that  the  theoretical  effects  of  the  limited  aperture  were  examined. 
The  effect  of  varying  the  crystal  angle  with  a  fixed  detector  angle  (with  respect  to  the 
electron  beam)  was  examined  The  variations  with  angle  of  I(n,0)/I{n,22.6®)  are  consistent 
with  theory  when  the  data  and  theory  are  normalized  to  a  given  angle  (in  this  case  22.6 
degrees )  However,  the  absolute  yield  still  greatly  exceeds  the  theoretical  values 
That  is,  for  both  the  far  and  the  near  regions  I(l,0)^*fM(l,0)theory»  \  i(n,0)/I(l,0)  is 
not  constant  with  respect  to  angle  as  is  the  theory.  The  maximum  in  this  ratio  occurs  at  9 
3.  If  current  is  not  monitored,  look  for  the  maximum  in  this  ratio  to  insure  the  crystal 
angle  is  63.  If  current  is  monitored  look  for  the  maximum  intensity  in  the  first  order  peak, 
1(1),  to  insure  the  crystal  angle  is  63. 

Additionally,  the  PXR  spectral  resolution  was  calculated  for  the  far  and  near  region. 
The  detector  subtended  a  solid  angle  of  2x10-4  sr.  for  the  far  region  and  2.37xlf^-3  sr  for 
the  near  region.  For  the  far  region  the  PXR  spectral  resolution  was  in  agreement  with 
theoiy  for  the  lower  orders,  but  was  narrower  than  expected  for  the  higher  orders  For 
the  near  region  the  PXR  spectral  resolution  was  in  agreement  with  theory,  such  that  the 
resolution  was  wider  for  a  given  energy  in  the  near  region  than  in  the  far  region. 
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Further  studies  need  to  be  done  to  fiirthe;  examine  the  effects  of  multiple  scattering  and 
mosaic  spread  By  comparing  a  thick  pure  crystal  such  as  Si  with  other  mosaic  crystals 
such  as  carbon,  the  effect  of  mosaic  spread  can  be  clarified  further. 

This  study  has  also  shown  that  although  alignment  is  not  extremely  critical  it  is 
important  to  take  great  care  in  focus  and  steering  of  the  electron  beam.  Additionally,  the 
effects  of  detector  drift,  or  shifts  due  to  other  causes  should  be  carefully  examined.  Many 
short  data  runs  sliould  be  collected  These  can  then  be  added  on  a  spreadsheet  after 
shifting  the  known  calibration  peaks  centroids  so  that  they  overlap,  thus  minimizing  these 
effects. 

It  is  extremely  important  for  calibration  purposes  to  have  an  x-ray  fluorescent  foil 
behind  the  crystal.  This  not  only  greatly  improved  the  analysis,  it  allowed  for  the 
calculation  of  the  LINAC  current  and  the  subsequent  calculation  of  the  absolute  number 
of  photons/electron.  Care  should  be  taken  in  selecting  a  fluorescent  foil  whose  x-ray  de¬ 
excitation  energy  falls  between  two  PXR  spectral  peaks.  Tin  worked  fairly  well  for 
carbon  at  shallower  angles,  but  overlapped  the  5th  order  peak  near  the  Bragg  angle, 
complicating  the  analysis. 

PXR  from  a  thick  mosaic  carbon  graphite  crystal  has  been  shown  to  be  a  spectrally 
bright  source  of  narrow  bandwidth  at  energies  equal  to  approximately  5  KeV  and  integer 
multiples  thereof  Further  study  with  other  crystals  could  provide  additional  sources  of 
PXR  at  other  energies  as  well. 


86 


APPENDIX  A.  FIGURES  OF  ROCKING  DATA 


Figure  42.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  5  June 
at  an  angle  of  19, 1  degrees.  This  spectrum  shows  negligible  PXR  peaks  and  a  tin  peak  at 
about  25  keV. 
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Figure  47.  PXR  spectrum  from  the  {002}  planes  of  carbon  gr^hite,  obtained  on  5  June 
at  an  angle  of  22,6  degrees.  This  spectrum  shows  six  orders  of  PXR  peaks.  The  fifth 
order  peak  and  the  Tin  calibration  peak  are  overlapping,  with  the  Tin  peak  at  a  higher 
energy  than  the  n=5  peak. 


Figure  48.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  5  June 
at  an  angle  of  23.1  degrees.  This  spectrum  shows  five  orders  of  PXR  peaks.  The  fifth 
order  peak  and  the  Tin  calibration  peak  are  overlapping  even  more  than  at  22.6  degrees, 
with  the  Tin  peak  at  a  higher  energy  than  the  n=5  peak. 
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Figure  S3.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  20. 1  degrees.  This  spectrum  shows  two  orders  of  PXR  peaks  and  the  Tin 
K^andKj)  peaks 
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Figure  54.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  21.1  degrees.  This  spectrum  shows  the  first  and  second  order  PXR  peak 
and  the  Tin  and  peaks. 
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Figure  55.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  21.6  degrees.  This  spectrum  shows  the  first  3  orders  of  PXR  peaks  and  the 
Tin  and  Kp  peaks. 
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Figure  56.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  22. 1  degrees.  This  spectrum  shows  four  orders  of  PXR  peaks  and  the  Tin 
and  Kp  peaks. 


101 


350i 


Figure  57.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  22.6  degrees.  This  spectrum  shows  five  orders  of  PXR  peaks  with  the  Tin 
peak  and  n=5  peaks  overlapping. 
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Figure  58.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtain^  on  22  July 
at  an  angle  of  23.1  degrees.  This  spectrum  shows  five  orders  of  FXR  peaks  The  fifth 
order  peak  and  the  Tin  calibration  peak  are  overlapping. 
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Figure  59.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  23. 6  degrees.  This  spearum  shows  possible  five  orders  of  PXR  peaks  The 
fifth  order  peak  and  the  Tin  calibration  peak  are  overlapping  . 
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Figure  60.  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  24  1  degrees.  This  spectrum  shows  possibly  four  orders  of  PXR  peaks  and 
a  tin  peak  at  about  25  keV. 
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Figure  61 .  PXR  spectrum  from  the  {002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  24,6  degrees.  This  spectrum  shows  three  orders  of  PXR  peaks  and  a  tin 
peak  at  about  25  keV 
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Figure  62.  PXR  spectrum  from  the  (002}  planes  of  carbon  graphite,  obtained  on  22  July 
at  an  angle  of  25.6  degrees.  This  spectrum  shows  two  orders  of  PXR  peaks  and  the  Tin 
and  Kp  peaks. 
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APPENDIX  B 


A.  PROGRAM  INC 

This  program  was  written  in  C  to  convert  the  ASCII  output  from  the  PHA-II  software 
to  a  file  readable  by  QPRO,  It  binned  every  four  channels  into  one  bin  and  assigned  the 
value  of  the  third  bin's  energy  to  that  bin.  This  program  allowed  the  8192  channels  to  be 
displayed  in  2048  channels  for  comparison  with  the  spectra  obtained  using  the  2048 
channel  PHA-I  software.  To  prepare  the  file,  an  editor  was  used  to  eliminate  the  commas, 
remove  the  heading  comments  and  add  "2045"  to  the  top  of  the  file  The  output  was  in 
the  form  of  "New  channel  number.  Energy,  counts"  (comma's  included)  The  parameter 
for  the  conversion  from  channel  to  energy  must  be  changed  from  experiment  to 
experiment  depending  upon  the  calibration  for  that  data.  See  Chapter  III  section  C 
equations  (28)-(30). 

#include  <stdio.h> 

#include  <math.h> 

main()  {int  ea,  eb,  ec,  ed,  i,  N,  counts; 
float  ca,  cb,  cc,  cd,  E; 
int  roia,  roib,  roic,  roid; 
scanf("%d",&n); 

E=0.0; 

counts=0.0;  ca=0;  cb=0;  cc*0;  cd=0; 

for  (i=0;i<N;++i)  {  scanfi["%d%f%d",&ea,&ca,&roia); 

scanf("%d%P/od",&cb,&cb,&roib); 

scanfl;"%d%f%d",&ec,&cc,&roic); 

scan^"%d%P/ttd",&ed,&cd,&roid); 

counts*(intXca+cb+cc+cd); 
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E=0.00472*(i*4+2)-1.04; 
printfC'%d'’,i);  printf(",''); 
printf(”%g",E);  printfl","), 
printf("%d'’.counts);  printf("\n'');  }  } 

B.  COMPUTER  INTEGRATION  PROGRAM  EXAMPLE 

This  program  was  used  to  determine  the  absolute  intensity  and  the  energy  bandwidth 
for  the  PXR  spectra.  The  various  parameters  changed  based  upon  the  experimental 
conditions.  This  example  uses  the  variable  for  the  far  region  at  a  crystal  angle  of  22  5 
degrees  It  includes  the  effects  of  multiple  scattering  and  mosaic  spread 
n=l 

enbeam=90.0 

nvec=(2.01, 1 .34,0.0946,0  624,0.402,0.268) 

fl=flvec[n] 

z=6 

an=12 

thetamm=0  4 
thetabb=22.5 
dang=6.708 

Uvec=  (0.02 1 3,0. 1 8 1 ,0. 573, 1 . 1 04, 1 . 586, 1 . 939) 

La=Lavec[n] 

dthetax=l5/l002.l9 

dthetay*=15/1002.19 

dnm=dang/10 

gam=(.  5 1 1  +enbeam)/.  5 1 1 

thetam*(pi/ 1 80)*thetamm 

thctab*(pj/ 1 80)*thetabb 
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L=  139/sin(thetab) 

num=La'^2-(l  *  !a+  La''3  )•  exp(  -L/La) 

den=La*(  1  -exp(-L/La) 

Tave=nuni/den 

k=0. 1 57*z*(z+l  )*Tave/An*enb€ain''2) 

Thetassqr“k*Ln(  1 . 13*z'^(4/3)an'^-l  )*Tave) 

energy=n*2*pi*(197.32)/(dnm*sin(thet^)) 

lambcm=0. 000 1 24/energy 

lambang=lambcm*  1 00000000 

x0=78.192*6/(encrgy)^2 

X 1 0sqr={  1  039e- 1 2)*(lambang^4)*fl ''2 

thetap=sqrt((  l/gam''2)+xo-*-thetassqr->-thetam^2) 

alphay=dthctay/(2*thetap) 

alphax=dthetax/(2*thctap) 

cs2=(cos(2  *thetab))''2 

j  1  (x)=(  1  +x^2  *(  1  +cs2))/(  1  +x''2r  1 . 5 

Q(x)*atan(alphay/(  1  +x''2y'.  5) 

j2  3(x)=x^2*cs2*((aIphay/(  I  +x''2))*  1  /(alphay^2+ 1  +x''2)+Q(x)/(  1  +x''2r  15) 
j2(x)=(j2a(x)+Q(x)/(  1  +x^2r  5.alphay/(alphay^2+ 1  +x^2))) 
factor=(  1  -cxp(-L/La)) 

coefF=xl0sqr*La/(4*sin(thetab)*cos(thctab)!ainbcm*thetap*energy*  1 37)*factor 

dNd£  1  (x>=pi  *cocff*j  1  (x) 

dNdE2(x)«2*coeff*j2(x) 

xdNdE2(x)*x*dNdE2(x) 

xH(pi/l  80)*(45-2*thetabb)+.  5dthetax)/thetap 

x2“((p«/l  80)*(45-2*thetaW)>-.5dthetax)/tIictap 

xccnt«»(iii({le-5*xdNdE2(x)),x*xl  to  x2))/(in((le-5*dNdE2(x)),x*xl  to  x2)) 


no 


ecent=xcent*energy*thetap/tan(thetab)+energy 
area=in(dNdE2(x),  x  =  xl  to  x2) 
resolution  -(x  1  -x2)*energy*thetap/tan(thetab). 
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